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Abstract 

The development and applications of the method of quasiaverages developed by N. N. 
Bogoliubov to quantum statistical physics and to quantum solid state theory and, in particular, 
to quantum theory of magnetism, were analyzed. The problem of finding the ferromagnetic, 
antiferromagnetic and superconducting symmetry broken solutions of the correlated lattice 
fermion models was discussed within the irreducible Green functions method. A unified scheme 
for the construction of generalized mean fields (elastic scattering corrections) and self-energy 
(inelastic scattering) in terms of the Dyson equation was generalized in order to include the 
source fields. The interrelation of the Bogoliubov's idea of quasiaverages and the concepts of 
symmetry breaking and quantum protectorate was discussed briefly in the context of quantum 
statistical physics. The idea of quantum protectorate reveals the essential difference in the 
behaviour of the complex many-body systems at the low-energy and high-energy scales. It was 
shown that the role of symmetry (and the breaking of symmetries) in combination with the 
degeneracy of the system was reanalyzed and essentially clarified within the framework of the 
method of quasiaverages. The complementary notion of quantum protectorate might provide 
distinctive signatures and good criteria for a hierarchy of energy scales and the appropriate 
emergent behavior. 
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1 Introduction 



It is well known that symmetry principles play a crucial role in physicsP^ The theory of symmetry 
is a basic tool for understanding and formulating the fundamental notions of physics P According 
to F. WilczekS "the primary goal of fundamental physics is to discover profound concepts that 
illuminate our understanding of nature" . It is known that symmetry is a driving force in the shap- 
ing of physical theory; moreover, the primary lesson of physics of last century is that the secret 
of nature is symmetry. Every symmetry leads to a conservation law; the well known examples 
are the conservation of energy, momentum and electrical charge. A variety of other conservation 
laws can be deduced from symmetry or invariance properties of the corresponding Lagrangian or 
Hamiltonian of the system. According to Noether theorem, every continuous symmetry transfor- 
mation under which the Lagrangian of a given system remains invariant implies the existence of 
a conserved function. 2 Many fundamental laws of physics in addition to their detailed features 
possess various symmetry properties. These symmetry properties lead to certain constraints and 
regularities on the possible properties of matter. Thus the principles of symmetries belong to the 
underlying principles of physics. 

It is known that when the Hamiltonian of a system is invariant under a symmetry operation, but 
the ground state is not, the symmetry of the system can be spontaneously brokenP Symmetry 
breaking is termed spontaneous when there is no explicit term in a Lagrangian which manifestly 
breaks the symmetry. Symmetries and breaking of symmetries play an important role in statistical 
physics, quantum field theory, physics of elementary particles, etcP^ 

In physics, spontaneous symmetry breaking occurs when a system that is symmetric with respect 
to some symmetry group goes into a vacuum state that is not symmetric. When that happens, the 
system no longer appears to behave in a symmetric manner. It is a phenomenon that naturally 
occurs in many situations. The mechanism of spontaneous symmetry breaking is usually under- 
stood as the mechanism responsible for the occurrence of asymmetric states in quantum systems 
in the thermodynamic limit and is used in various field of quantum physicsP The intriguing 
mechanism of spontaneous symmetry breaking is a unifying concept that lie at the basis of most 
of the recent developments in theoretical physics, from statistical mechanics to many-body theory 
and to elementary particles theoryP^ 

It should be stressed that symmetry implies degeneracy. The greater the symmetry, the greater 
the degeneracy. The study of the degeneracy of the energy levels plays a very important role in 
quantum physics. It is of importance to emphasize that when spontaneous symmetry breaking 
takes place, the ground state of the system is degenerate. Substantial progress in the under- 
standing of the broken symmetry concept was connected with Bogoliubov's fundamental ideas on 
quasiaverages Studies of degenerate systems led Bogoliubov in 1960-61 to the formulation 
of the method of quasiaverages. This method has proved to be a universal tool for systems 
whose ground states become unstable under small perturbations. Thus the role of symmetry (and 
the breaking of symmetries) in combination with the degeneracy of the system was reanalyzed 
and essentially clarified by N. N. Bogoliubov in 1960-1961. He invented and formulated a powerful 
innovative idea of quasiaverages in statistical mechanics JSEMffllMISllIS The very elegant work of 
N. N. Bogoliubov^ has been of great importance for a deeper understanding of phase transitions, 
superfluidity and superconductivity, quantum theory of magnetism^ and other fields of equilib- 
rium and nonequilibrium statistical rnechanics.^^HIllElIIElliQHlS] ^he concept of quasiaverages is 
indirectly related to the theory of phase transition. The instability of thermodynamic averages 
with respect to perturbations of the Hamiltonian by a breaking of the invariance with respect 
to a certain group of transformations means that in the system transition to an extremal state 
occurs. The mathematical apparatus of the method of quasiaverages includes the Bogoliubov 
theorerr F^^ ^^ on singularities of type 1/g 2 and the Bogoliubov inequality for Green and cor- 
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relation functions as a direct consequence of the method. It includes algorithms for establishing 
non-trivial estimates for equilibrium quasiaverages, enabling one to study the problem of ordering 
in statistical systems and to elucidate the structure of the energy spectrum of the underlying 
excited states. Thus the Bogoliubov's idea of quasiaverages is an essential conceptual advance of 
modern physics. 

Many-particle systems where the interaction is strong have often complicated behavior, and re- 
quire nonperturbative approaches to treat their properties. Such situations are often arise in 
condensed matter systems. Electrical, magnetic and mechanical properties of materials are emer- 
gent collective behaviors of the underlying quantum mechanics of their electrons and constituent 
atoms. A principal aim of solid state physics and materials science is to elucidate this emergence. 
A full achievement of this goal would imply the ability to engineer a material that is optimum 
for any particular application. The current understanding of electrons in solids uses simplified 
but workable picture known as the Fermi liquid theory. This theory explains why electrons in 
solids can often be described in a simplified manner which appears to ignore the large repulsive 
forces that electrons are known to exert on one another. There is a growing appreciation that this 
theory probably fails for entire classes of possibly useful materials and there is the suspicion that 
the failure has to do with unresolved competition between different possible emergent behaviors. 
It is appropriate to note here that the emergent properties of matter were analyzed and discussed 
by R. Laughlin and D. Pines^EID f r0 m a general point of view (see also RefP^J). They introduced 
a unifying idea of quantum protectorate. This concept belongs also to the underlying principles 
of physics. The idea of quantum protectorate reveals the essential difference in the behavior of 
the complex many-body systems at the low-energy and high-energy scales. The existence of two 
scales, low-energy and high-energy, in the description of physical phenomena is used in physics, 
explicitly or implicitly. 

From the other hand, it was recognized for many years that the strong correlation in solids exist 
between the motions of various particles (electrons and ions, i.e. the fermion and boson degrees 
of freedom) which arise from the Coulomb forces. The most interesting objects are metals and 
their compounds. They are invariant under the translation group of a crystal lattice and have 
lattice vibrations as well as electron degrees of freedom. There are many evidences for the im- 
portance of many-body effects in these systems. Within the semi-phenomenological theory it was 
suggested that the low-lying excited states of an interacting Fermi gas can be described in terms 
of a set of "independent quasiparticles". However, this was a phenomenological approach and did 
not reveal the nature of relevant interactions. An alternative way of viewing quasiparticles, more 
general and consistent, is through the Green function scheme of many-body theory EMZ2112M22] n 
becomes clear that only a thorough experimental and theoretical investigation of quasiparticle 
many-body dynamics of the many particle systems can provide the answer on the relevant mi- 
croscopic picture.^ In our works, we discussed the microscopic view of a dynamic behaviour 
of various interacting many-body systems on a latticeP^SBS] A comprehensive description of 
transition and rare-earth metals and alloys and other materials (as well as efficient predictions of 
properties of new materials) is possible only in those cases, when there is an adequate quantum- 
statistical theory based on the information about the electron and crystalline structures. The 
main theoretical problem of this direction of research, which is the essence of the quantum theory 
of magnetism, is investigations and improvements of quantum-statistical models describing the 
behavior of the complex compounds and materials in order to take into account the main features 
of their electronic structure, namely, their dual "band-atomic" nature.^ The construction of a 
consistent theory explaining the electronic structure of these substances encounters serious diffi- 
culties when trying to describe the collectivization-localization duality in the behavior of electrons. 
This problem appears to be extremely important, since its solution gives us a key to understanding 
magnetic, electronic, and other properties of this diverse group of substances. The author of the 
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present paper investigated the suitability of the basic models with strong electron correlations and 
with a complex spectrum for an adequate and correct description of the dual character of electron 
statespZl A universal mathematical formalism was developed for this investigationPSEH n takes 
into account the main features of the electronic structure and allows one to describe the true 
quasiparticle spectrum, as well as the appearance of the magnetically ordered, superconducting, 
and dielectric (or semiconducting) states. With a few exceptions, diverse physical phenomena 
observed in compounds and alloys of transition and rare-earth metals, cannot be explained in 
the framework of the mean- field approximation, which overestimates the role of inter-electron 
correlations in computations of their static and dynamic characteristics. The circle of questions 
without a precise and definitive answer, so far, includes such extremely important (not only from 
a theoretical, but also from a practical point of view) problems as the adequate description of 
quasiparticle dynamics for quantum-statistical models in a wide range of their parameter values. 
The source of difficulties here lies not only in the complexity of calculations of certain dynamic 
properties (such as, the density of states, electrical conductivity, susceptibility, electron-phonon 
spectral function, the inelastic scattering cross section for slow neutrons), but also in the ab- 
sence of a well-developed method for a consistent quantum-statistical analysis of a many-particle 
interaction in such systems. A self-consistent field approach was used in the papers^E^SS] f or 
description of various dynamic characteristics of strongly correlated electronic systems. It allows 
one to consistently and quite compactly compute quasiparticle spectra for many-particle systems 
with strong interaction taking into account damping effects. The correlation effects and quasipar- 
ticle damping are the determining factors in analysis of the normal properties of high-temperature 
superconductors, heavy fermion compounds, etc. We also formulated a general scheme for a the- 
oretical description of electronic properties of many-particle systems taking into account strong 
inter-electron correlations! 22 ^ 38 ^ The scheme is a synthesis of the method of two-time tempera- 
ture Green's functionaries] anc i th e diagram technique. An important feature of this approach is 
a clear-cut separation of the elastic and inelastic scattering processes in many-particle systems 
(which is a highly nontrivial task for strongly correlated systems). As a result, one can construct a 
correct basic approximation in terms of generalized mean fields (the elastic scattering corrections), 
which allows one to describe magnetically ordered or superconducting states of the system. The 
residual correlation effects, which are the source of quasiparticle damping, are described in terms 
of the Dyson equation with a formally exact representation for the mass operator. 
In the present paper we will discuss some applications of the symmetry principles to quantum 
and statistical physics and quantum solid state theory in the light of our results on quasiparticle 
many-body dynamics. 

2 Bogoliubov's Quasiaverages in Statistical Mechanics 

In the work of N. N. Bogoliubov "Quasiaverages in Problems of Statistical Mechanics" the innova- 
tive notion of quasiavereg^^ was introduced and applied to various problem of statistical physics. 
In particular, quasiaverages of Green's functions constructed from ordinary averages, degeneration 
of statistical equilibrium states, principle of weakened correlations, and particle pair states were 
considered. In this framework the l/t/ 2 -type properties in the theory of the superfluidity of Bose 
and Fermi systems, the properties of basic Green functions for a Bose system in the presence of 
condensate, and a model with separated condensate were analyzed. 

The method of quasiaverages is a constructive workable scheme for studying systems with spon- 
taneous symmetry breakdown. A quasiaverage is a thermodynamic (in statistical mechanics) or 
vacuum (in quantum field theory) average of dynamical quantities in a specially modified averag- 
ing procedure, enabling one to take into account the effects of the influence of state degeneracy of 
the system. The method gives the so-called macro-objectivation of the degeneracy in the domain 
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of quantum statistical mechanics and in quantum physics. In statistical mechanics, under sponta- 
neous symmetry breakdown one can, by using the method of quasiaverages, describe macroscopic 
observable within the framework of the microscopic approach. 

In considering problems of findings the eigenfunctions in quantum mechanics it is well known 
that the theory of perturbations should be modified substantially for the degenerate systems. In 
the problems of statistical mechanics we have always the degenerate case due to existence of the 
additive conservation laws. The traditional approach to quantum statistical mechanics^ is based 
on the unique canonical quantization of classical Hamiltonians for systems with finitely many 
degrees of freedom together with the ensemble averaging in terms of traces involving a statistical 
operator p. For an operator A corresponding to some physical quantity A the average value of A 
will be given as 

(A) H = TrpA; p = exp"^ /Trexp~ PH , (2-1) 

where H is the Hamiltonian of the system, j3 = 1/kT is the reciprocal of the temperature. 
The core of the problem lies in establishing the existence of a thermodynamic limit (such as 
N/V = const, V — > oo, N = number of degrees of freedom, V = volume) and its evaluation 
for the quantities of interest. Thus in the statistical mechanics the average (A) of any dynamical 
quantity A is defined in a single- valued way. In the situations with degeneracy the specific problems 
appear. In quantum mechanics, if two linearly independent state vectors (wavefunctions in the 
Schroedinger picture) have the same energy, there is a degeneracy. In this case more than one 
independent state of the system corresponds to a single energy level. If the statistical equilibrium 
state of the system possesses lower symmetry than the Hamiltonian of the system (i.e. the situation 
with the symmetry breakdown), then it is necessary to supplement the averaging procedure ()2. 1|) 
by a rule forbidding irrelevant averaging over the values of macroscopic quantities considered 
for which a change is not accompanied by a change in energy. This is achieved by introducing 
quasiaverages, that is, averages over the Hamiltonian H v g supplemented by infinitesimally-small 
terms that violate the additive conservations laws H v g = H + u{e ■ M), (y — >■ 0). Thermodynamic 
averaging may turn out to be unstable with respect to such a change of the original Hamiltonian, 
which is another indication of degeneracy of the equilibrium state. According to Bogoliubovp^ 
the quasiaverage of a dynamical quantity A for the system with the Hamiltonian H u g is defined 
as the limit 

< A^=lim(A) u g, (2.2) 

v — ^0 

where {A) u g denotes the ordinary average taken over the Hamiltonian H u g, containing the small 
symmetry-breaking terms introduced by the inclusion parameter v, which vanish as v — > after 
passage to the thermodynamic limit V — > oo. It is important to note that in this equation limits 
cannot be interchanged. Thus the existence of degeneracy is reflected directly in the quasiaverages 
by their dependence upon the arbitrary unit vector e. It is also clear that 

(A) = J <A> de. (2.3) 

According to definition (12.3p . the ordinary thermodynamic average is obtained by extra averaging 
of the quasiaverage over the symmetry-breaking group. Thus to describe the case of a degenerate 
state of statistical equilibrium quasiaverages are more convenient, more physical, than ordinary 
averages.^ The latter are the same quasiaverages only averaged over all the directions e. 
It is necessary to stress, that the starting point for Bogoliubov's workP^ was an investigation of 
additive conservation laws and selection rules, continuing and developing the approach by P. Curie 
for derivation of selection rules for physical effects (see also RefP). Bogoliubov demonstrated that 
in the cases when the state of statistical equilibrium is degenerate, as in the case of the Heisenberg 
ferromagnet, one can remove the degeneracy of equilibrium states with respect to the group of 
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spin rotations by including in the Hamiltonian H an additional noninvariant term vM z V with an 
infinitely small v. For the Heisenberg ferromagnet the ordinary averages must be invariant with 
regard to the spin rotation group. The corresponding quasiaverages possess only the property of 
covariance. Thus the quasiaverages do not follow the same selection rules as those which govern 
ordinary averages, due to their invariance with regard to the spin rotation group. It is clear that 
that the unit vector e, i.e., the direction of the magnetization M vector, characterizes the degen- 
eracy of the considered state of statistical equilibrium. In order to remove the degeneracy one 
should fix the direction of the unit vector e. It can be chosen to be along the z direction. Then 
all the quasiaverages will be the definite numbers. This is the kind that one usually deals with in 
the theory of ferromagnetism. 

The value of the quasi-average (12. 2D may depend on the concrete structure of the additional term 
AH = H u — H, if the dynamical quantity to be averaged is not invariant with respect to the 
symmetry group of the original Hamiltonian H. For a degenerate state the limit of ordinary 
averages (|2.3p as the inclusion parameters v of the sources tend to zero in an arbitrary fashion, 
may not exist. For a complete definition of quasiaverages it is necessary to indicate the manner 
in which these parameters tend to zero in order to ensure convergence.^ On the other hand, in 
order to remove degeneracy it suffices, in the construction of H, to violate only those additive 
conservation laws whose switching lead to instability of the ordinary average. Thus in terms of 
quasiaverages the selection rules for the correlation functions^^ that are not relevant are those 
that are restricted by these conservation laws. 

By using H u , we define the state uj{A) = (A) u and then let v tend to zero (after passing to the 
thermodynamic limit) PEUE] if a ll 

averages uj (A) get infinitely small increments under infinitely 
small perturbations v, this means that the state of statistical equilibrium under consideration is 
nondegeneratePHE' However, if some states have finite increments as v — > 0, then the state is 
degenerate. In this case, instead of ordinary averages {A)h, one should introduce the quasiaver- 
ages (I2.2h . for which the usual selection rules do not hold. 

The method of quasiaverages is directly related to the principle weakening of the correlatiorfl^^^ 
in many-particle systems. According to this principle, the notion of the weakening of the corre- 
lation, known in statistical mechanics pSHHEHl i n the- case of state degeneracy must be interpreted 
in the sense of the quasiaverages P 

The quasiaverages may be obtained from the ordinary averages by using the cluster property 
which was formulated by Bogoliubov.H' This was first done when deriving the Boltzmann equa- 
tions from the chain of equations for distribution functions, and in the investigation of the model 
Hamiltonian in the theory of superconductivityP^^^^ To demonstrate this let us consider 
averages (quasiaverages) of the form 

F(ti,x 1) ...t n ,x n ) = {...&(t 1 ,xi)...-$(tj,x j )...), (2.4) 

where the number of creation operators ^ may be not equal to the number of annihilation opera- 
tors We fix times and split the arguments (t\, x\, . . . t n ,x n ) into several clusters (. . . , t a ,x a , ...),..., 
(. . . ,tp,xp, . . .). Then it is reasonably to assume that the distances between all clusters \x a — xp\ 
tend to infinity. Then, according to the cluster property, the average value (j2.4|) tends to the prod- 
uct of averages of collections of operators with the arguments (. . . , t a ,x a , ...),..., (. . . , tp,xp, . . .) 

lim F(ti,xi,...t n ,x n ) = F(...,t a ,x a ,...)...F(...,t/3,Xi3,...). (2.5) 

Xg |— too 

For equilibrium states with small densities and short-range potential, the validity of this property 
can be proved.^ For the general case, the validity of the cluster property has not yet been proved. 
Bogoliubov formulated it not only for ordinary averages but also for quasiaverages, i.e., for anoma- 
lous averages, too. It works for many important models, including the models of superfluidity^^ 
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and superconductivitjPESEEl ( see also RefsPHSl). 

To illustrate this statement consider Bogoliubov's theory of a Bose-system with separated con- 
densate, which is given by the Hamiltoniarl^"^ 

H A = [ yU x )(- — )V(x)dx - ii ! &(x)V(x)dx (2.6) 
Ja 2m Ja 

+7; I ^{ x i)^H x 2)<£>(xi ~ x 2 )^{x 2 )^{xi)dxidx 2 . 

2 JA2 

This Hamiltonian can be written also in the following form 

H A = H + H 1= I y\ q )(-A)y( q )dq (2.7) 



+ - / y\q)&(q')<5>(q-q')^(q')y(q)dqdq'. 
' I A 2 
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Here, ^(q), and W(q) are the operators of annihilation and creation of bosons. They satisfy the 
canonical commutation relations 

[*(<?), = Kq - <?'); *(</)] = *V)] = 0. (2.8) 

The system of bosons is contained in the cube A with the edge L and volume V. It was assumed 
that it satisfies periodic boundary conditions and the potential $((?) is spherically symmetric and 
proportional to the small parameter. It was also assumed that, at temperature zero, a certain 
macroscopic number of particles having a nonzero density is situated in the state with momentum 
zero. The operators ^(q), and ^(q) can be represented in the form 

= a /W; &(q) = 4/W, (2.9) 

where ciq and are the operators of annihilation and creation of particles with momentum zero. 
To explain the phenomenon of superfluidity, one should calculate the spectrum of the Hamilto- 
nian, which is quite a difficult problem. Bogoliubov suggested the idea of approximate calculation 
of the spectrum of the ground state and its elementary excitations based on the physical nature 
of superfluidity. His idea consists of a few assumptions. The main assumption is that at temper- 
ature zero the macroscopic number of particles (with nonzero density) has the momentum zero. 
Therefore, in the thermodynamic limit, the operators ao/W and a^/VV commute 

lim 

y-»oo 



a /W,al/W]=^^0 (2.10) 



and are c-numbers. Hence, the operator of the number of particles No = a^ao is a c-number, 
too. It is worth noting that the Hamiltonian (|2.7h is invariant under the gauge transformation 
dk = exp(iip)dk, at = exp(— i<p)a\, where cp is an arbitrary real number. Therefore, the averages 
(ao/y/V) and (o,q/\/V) must vanish. But this contradicts to the assumption that clq/VV and 
a /\/V must become c-numbers in the thermodynamic limit. In addition it must be taken into 
account that a^ao/V = Nq/V ^ which implies th&tao/VV = No ex.p(ia)/W ^ and a\/y/V = 
No exp(— ia)/\/V ^ 0, where a is an arbitrary real number. This contradiction may be overcome 
if we assume that the eigenstates of the Hamiltonian are degenerate and not invariant under gauge 
transformations, i.e., that a spontaneous breaking of symmetry takes place. 

Thus the averages {ao/y/V) and (aJ/VF), which are nonzero under spontaneously broken gauge 
invariance, are called anomalous averages or quasiaverages. This innovative idea of Bogoliubov 
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penetrate deeply into the modern quantum physics. The systems with spontaneously broken 
symmetry are studied by use of the transformation of the operators of the form 

= a /W + 8(qy, &(q) = al/W + 0*(q), (2.11) 

where ao/VV and a\/\/V are the numbers first introduced by Bogoliubov in 1947 in his inves- 
tigation of the phenomenon of superfluidityP^^^ The main conclusion was made that for the 
systems with spontaneously broken symmetry, the quasiaverages should be studied instead of 
the ordinary averages. It turns out that the long-range order appears not only in the system of 
Bose-particles but also in all systems with spontaneously broken symmetry. Bogoliubov's papers 
outlined above anticipated the methods of investigation of systems with spontaneously broken 
symmetry for many years. 

As mentioned above, in order to explain the phenomenon of superfluidity, Bogoliubov assumed 
that the operators ao/VV and a^/VV become c-numbers in the thermodynamic limit. This state- 
ment was rigorously proved in the papers by Bogolyubov and some other authors. Bogolyubov's 
proof was based on the study of the equations for two-time Green's functions^ and on the as- 
sumption that the cluster property holds. It was proved that the solutions of equations for Green's 
functions for the system with Hamiltonian (|2.7p coincide with the solutions of the equations for 
the system with the same Hamiltonian in which the operators ao/VV and a /W are replaced by 
numbers. These numbers should be determined from the condition of minimum for free energy. 
Since all the averages in both systems coincide, their free energies coincide, too. 
It is worth noting that the validity of the replacement of the operators ao and aj by c-numbers 
in the thermodynamic limit was confirmed in the numerous subsequent publications of various 
authors. Thus Bogoliubov's 1947 analysis of the many-body Hamiltonian by means of a c-number 
substitution for the most relevant operators in the problem, the zero-momentum mode operators, 
was justified rigorously. Since the Bogoliubov's 1947 analysis is one of the key developments in 
the theory of the Bose gas, especially the theory of the low density gases currently at the forefront 
of experiment H3HS1 this result is of importance for the legitimation of that theory. Additional 
arguments were given in study, where the Bose-Einstein condensation and spontaneous £7(1) sym- 
metry breaking were investigated on the basis of Bogoliubov's truncated Hamiltonian Hb for a 
weakly interacting Bose system, and adding a £7(1) symmetry breaking term y/V(\ao + A*Oq) 
to Hb, It was shown also, by using the coherent state theory and the mean- field approximation 
rather than the c-number approximations, that the Bose-Einstein condensation occurs if and only 
if the £7(1) symmetry of the system is spontaneously broken. The real ground state energy and 
the justification of the Bogoliubov c-number substitution were given by solving the Schroedinger 
eigenvalue equation and using the self-consistent condition. Thus the Bogoliubov c-number sub- 
stitutions were fully correct and the symmetry breaking causes the displacement of the condensate 
state. 

The concept of quasiaverages was introduced by Bogoliubov on the basis of an analysis of many- 
particle systems with a degenerate statistical equilibrium state. Such states are inherent to various 
physical many-particle systems Those are liquid helium in the superfluid phase, metals in the 
superconducting state, magnets in the ferromagnetically ordered state, liquid crystal states, the 
states of superfluid nuclear matter, etc. (for a review, see In case of superconductivity, 

the source v ^ fc v(k)(a\^a_ k < + a_^a^) was inserted in the BCS-Bogoliubov Hamiltonian, and 
the quasiaverages were defined by use of the Hamiltonian H v . In the general case, the sources are 
introduced to remove degeneracy. If infinitesimal sources give infinitely small contributions to the 
averages, then this means that the corresponding degeneracy is absent, and there is no reason to 
insert sources in the Hamiltonian. Otherwise, the degeneracy takes place, and it is removed by 
the sources. The ordinary averages can be obtained from quasiaverages by averaging with respect 
to the parameters that characterize the degeneracy. 
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N. N. Bogoliubov, Jr. 15 considered some features of quasiaverages for model systems with four- 
fermion interaction. He discussed the treatment of certain three-dimensional model systems which 
can be solved exactly. For this aim a new effective way of defining quasiaverages for the systems 
under consideration was proposed. 

Peletminskii and SokolovskiP^ have found general expressions for the operators of the flux den- 
sities of physical variables in terms of the density operators of these variables. The method of 
quasiaverages and the expressions found for the flux operators were used to obtain the averages 
of these operators in terms of the thermodynamic potential in a state of statistical equilibrium of 
a superfluid liquid. 

Vozyakov 4 ^ reformulated the theory of quantum crystals in terms of quasiaverages. He analyzed 
a Bose system with periodic distribution of particles which simulates an ensemble in which the 
particles cannot be regarded as vibrating independently about a position of equilibrium lattice 
sites. With allowance for macroscopic filling of the states corresponding to the distinguished sym- 
metry, a calculation was made of an excitation spectrum in which there exists a collective branch 
of gapless type. 

Peregoudov 4 ^ discussed the effective potential method, used in quantum field theory to study 
spontaneous symmetry breakdown, from the point of view of Bogoliubov's quasiaveraging proce- 
dure. It was shown that the effective potential method is a disguised type of this procedure. The 
catastrophe theory approach to the study of phase transitions was discussed and the existence 
of the potentials used in that approach was proved from the statistical point of view. It was 
shown that in the ease of broken symmetry, the nonconvex effective potential is not a Legendre 
transform of the generating functional for connected Green's functions. Instead, it is a part of 
the potential used in catastrophe theory. The relationship between the effective potential and 
the Legendre transform of the generating functional for connected Green's functions is given by 
Maxwell's rule. A rigorous rule for evaluating quasiaveraged quantities within the framework of 
the effective potential method was established. 

N. N. Bogoliubov, Jr. with M. Yu. Kovalevsky and co-authors^ developed a statistical approach 
for solving the problem of classification of equilibrium states in condensed media with sponta- 
neously broken symmetry based on the quasiaverage concept. Classification of equilibrium states 
of condensed media with spontaneously broken symmetry was carried out. The generators of 
residual and spatial symmetries were introduced and equations of classification for the order pa- 
rameter has been found. Conditions of residual symmetry and spatial symmetry were formulated. 
The connection between these symmetry conditions and equilibrium states of various media with 
tensor order parameter was found out. An analytical solution of the problem of classification of 
equilibrium states for superfluid media, liquid crystals and magnets with tensor order parameters 
was obtained. Superfluid 3 He, liquid crystals, quadrupolar magnetics were considered in detail. 
Possible homogeneous and heterogeneous states were found out. Discrete and continuous ther- 
modynamic parameters, which define an equilibrium state, allowable form of order parameter, 
residual symmetry, and spatial symmetry generators were established. This approach, which is 
alternative to the well-known Ginzburg-Landau method, does not contain any model assumptions 
concerning the form of the free energy as functional of the order parameter and does not employ 
the requirement of temperature closeness to the point of phase transition. For all investigated 
cases they found the structure of the order parameters and the explicit forms of generators of 
residual and spatial symmetries. Under the certain restrictions they established the form of the 
order parameters in case of spins 0, 1/2, 1 and proposed the physical interpretation of the studied 
degenerate states of condensed media. 
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3 Quantum Protectorate 



The "quantum protectorate" concept was formulated in paperP^l Its inventors, R. Laughlin and 
D. Pines, discussed the most fundamental principles of matter description in the widest sense of 
this word. They formulated their main thesis: emergent physical phenomena, which are regulated 
by higher physical principles, have a certain property, typical for these phenomena only. This 
property is their insensitivity to microscopic description. For instance, the crystalline state is 
the simplest known example of a quantum protectorate, a stable state of matter whose generic 
low-energy properties are determined by a higher organizing principle and nothing else. There are 
many other examples.^ These quantum protectorates, with their associated emergent behavior, 
provide us with explicit demonstrations that the underlying microscopic theory can easily have 
no measurable consequences whatsoever at low energies. The nature of the underlying the- 
ory is unknowable until one raises the energy scale sufficiently to escape protection. 
The existence of two scales, the low-energy and high-energy scales, relevant to the description 
of magnetic phenomena was stressed by the author of the present work in the papersj^E^l 
which were devoted to comparative analysis of localized and band models of quantum theory of 
magnetism. It was shown there, that the low-energy spectrum of magnetic excitations in the 
magnetically-ordered solid bodies corresponds to a hydrodynamic pole (k, cj — > 0) in the general- 
ized spin susceptibility x> which is present in the Heisenberg, Hubbard, and the combined s — d 
model. In the Stoner band model the hydrodynamic pole is absent, there are no spin waves there. 
At the same time, the Stoner single-particle's excitations are absent in the Heisenberg model's 
spectrum. The Hubbard model with narrow energy bands contains both types of excitations: the 
collective spin waves (the low-energy spectrum) and Stoner single-particle's excitations (the high- 
energy spectrum). This is a big advantage and flexibility of the Hubbard model in comparison 
to the Heisenberg model. The latter, nevertheless, is a very good approximation to the realistic 
behavior in the limit k,u —> 0, the domain where the hydrodynamic description is applicable, 
that is, for long wavelengths and low energies. The quantum protectorate concept was applied 
to the quantum theory of magnetism by the the present author in paperP where a criterion of 
applicability of models of the quantum theory of magnetism to description of concrete substances 
was formulated. The criterion is based on the analysis of the model's low-energy and high-energy 
spectra. There the detailed analysis was carried out of the idea of quantum protectorate^ in the 
context of quantum theory of magnetism.^ It was suggested that the difficulties in the formula- 
tion of quantum theory of magnetism at the microscopic level, that are related to the choice of 
relevant models, can be understood better in the light of the quantum protectorate concept.^ We 
argued that the difficulties in the formulation of adequate microscopic models of electron and mag- 
netic properties of materials are intimately related to dual, itinerant and localized behaviour of 
electrons.^ We formulated a criterion of what basic picture describes best this dual behaviour. 
The main suggestion was that quasiparticle excitation spectra might provide distinctive signatures 
and good criteria for the appropriate choice of the relevant model. A broad class of the problems 
of condensed matter physic d 22 ^ 23 ^ in the fields of the magnetism and superconductivity of complex 
materials were reconsidered in relation to these ideas. 



4 Irreducible Green Functions Method 

It was shown above that it becomes clear that only a thorough experimental and theoretical 
investigation of quasiparticle many-body dynamics of the many particle systems can provide the 
answer on the relevant microscopic picture. Many-particle systems where the interaction is strong 
have often complicated behavior, and require nonperturbative approaches to treat their properties. 
There are many different approaches to construction of generalized mean-field approximations; 
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however, all of them have a special-case character. The method of irreducible Green functions 
(IGF) allows one to tackle this problem in a more systematic fashion. In order to clarify this 
statement let us consider briefly the main ideas of the IGF approach that allows one to describe 
completely quasiparticle spectra with damping in a very natural way. When working with infinite 
hierarchies of equations for Green functions the main problem is finding the methods for their 
efficient decoupling, with the aim of obtaining a closed system of equations, which determine the 
Green functions. A decoupling approximation must be chosen individually for every particular 
problem, taking into account its character. This "individual approach" is the source of critique for 
being too not transparent, which sometimes appear in the papers using the causal Green functions 
and diagram technique. However, the ambiguities are also present in the diagram technique, when 
the choice of an appropriate approximation is made there. The decision, which diagrams one has 
to sum up, is obvious only for a narrow range of relatively simple problems. In the paper^ devoted 
to Bose-systems, and in the papers by the author of present work^EM^E] devoted to Fermi systems 
it was shown that for a wide range of problems in statistical mechanics and theory of condensed 
matter one can outline a fairly systematic recipe for constructing approximate solutions in the 
framework of irreducible Green's functions method. Within this approach one can look from a 
unified point of view at the main problems of fundamental characters arising in the method of 
two-time temperature Green functions. The method of irreducible Green functions is a useful 
reformulation of the ordinary Bogoliubov-Tyablikov method of equations of motionP^S 
We reformulated the two-time Green functions met ho j J^PoM UlMl ^ form which is especially 
adjusted to correlated fermion systems on a lattice and systems with complex spectra. A very 
important concept of the whole method is the generalized mean fields (GMFs), as it was formulated 
mJ22]|38] These GMFs have a complicated structure for a strongly correlated case and complex 
spectra, and are not reduced to the functional of mean densities of the electrons or spins when 
one calculates excitation spectra at finite temperatures. 

To clarify the foregoing, let us consider a retarded Green function of the formal 

G r = ((A(t),A\t'))) = -z8(t - t')([A(t)A^t%), V = ±1. (4.1) 

As an introduction to the concept of IGF, let us describe the main ideas of this approach in a 
symbolic and simplified form. To calculate the retarded Green function G(t — t'), let us write 
down the equation of motion for it 

uG(u) = ([A, A%) + (([A, | At)) w . (4.2) 

Here we use the notation ((A(t), A'(t'))) for the time-dependent Green function and ((A \ A^))^ 
for its Fourier transform.^ The notation [A, B]„ refers to commutation and anticommutation, 
depending on the value of r\ = ±. The essence of the method is as follows.^ It is based on the 
notion of the "IRRED U CIBLE" parts of Green functions (or the irreducible parts of the operators, 
A and A\ out of which the Green function is constructed) in terms of which it is possible, without 
recourse to a truncation of the hierarchy of equations for the Green functions, to write down the 
exact Dyson equation and to obtain an exact analytic representation for the self-energy operator. 
By definition, we introduce the irreducible part (ir) of the Green function 

W(([A,fl]_|At» = (([A,H]- - zA\A% (4.3) 

The unknown constant z is defined by the condition (or constraint) 

([[A,H]^\A%) = 0, (4.4) 

which is an analogue of the orthogonality condition in the Mori formalism EHEU us emphasize 
that due to the complete equivalence of the definition of the irreducible parts for the Green 
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functions { {ir) (([A, J5T]_|A+») and operators (^[A, = (L4, #]_) (ir) we will use both the 

notation freely ( is the same as <((A)( ir )|.B)) ). A choice one notation over another is 

determined by the brevity and clarity of notation only. From the condition (|4.4p one can find 

_ ([[A,H]^A%) _ My 

([A,A%) "Mo" (4 ' 5j 

Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, the 
irreducible Green functions are defined so that they cannot be reduced to the lower-order ones by 
any kind of decoupling. It is worth noting that the term "irreducible" in a group theory means a 
representation of a symmetry operation that cannot be expressed in terms of lower dimensional 
representations. Irreducible (or connected ) correlation functions are known in statistical mechan- 
ics. In the diagrammatic approach, the irreducible vertices are defined as graphs that do not 
contain inner parts connected by the G°-line. With the aid of the definition (14.3|) these concepts 
are expressed in terms of retarded and advanced Green functions. The procedure extracts all 
relevant (for the problem under consideration) mean-field contributions and puts them into the 
generalized mean-field Green function which is defined here as 

G0( W ) = y» At U ( 4. 6) 

To calculate the IGF (([A, B\-(t), A^t'))) in (|4.2p , we have to write the equation of motion for 
it after differentiation with respect to the second time variable t'. The condition of orthogonality 
(14. 4p removes the inhomogeneous term from this equation and is a very crucial point of the whole 
approach. If one introduces the irreducible part for the right-hand side operator, as discussed above 
for the "left" operator, the equation of motion (|4.2p can be exactly rewritten in the following form: 



G = G° + G°PG°. (4.7) 

The scattering operator P is given by 

P = (M )- 1 ( iir) (([A,HU[A\HU))^)(M )- 1 . (4.8) 

The structure of equation ( I4.8P enables us to determine the self-energy operator M by analogy 
with the diagram technique 

P = M + MG°P. (4.9) 

We used here the notation M for self-energy (mass operator in quantum field theory). From 
the definition (I4.9P it follows that the self-energy operator M is defined as a proper (in the 
diagrammatic language, "connected") part of the scattering operator M = (P) p - As a result, we 
obtain the exact Dyson equation for the thermodynamic double-time Green functions 

G = G° + G°MG. (4.10) 

The difference between P and M can be regarded as two different solutions of two integral equa- 
tions (|4,7p and (|4.10p . However, from the Dyson equation (|4.10p only the full GF is seen to be 
expressed formal solution of the form 

G= [(G )- 1 -M]-\ (4.11) 

Equation (|4.1ip can be regarded as an alternative form of the Dyson equation (|4.10p and the 
definition of M provides that the generalized mean-field GF G° is specified. On the contrary , for 
the scattering operator P, instead of the property + G M = 1, one has the property 

(G )- 1 -G- 1 = PG°G-\ 
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Thus, the very functional form of the formal solution (|4.1ip precisely determines the difference 
between P and M. 

Thus, by introducing irreducible parts of GF (or irreducible parts of the operators, out of which 
the GF is constructed) the equation of motion (|4.2p for the GF can exactly be (but using the 
orthogonality constraint (|4.4p ) transformed into the Dyson equation for the double-time thermal 
GF (I4.10p . This result is very remarkable because the traditional form of the GF method does 
not include this point. Notice that all quantities thus considered are exact. Approximations 
can be generated not by truncating the set of coupled equations of motions but by a specific 
approximation of the functional form of the mass operator M within a self-consistent scheme 
expressing M in terms of the initial GF 

M « F[G]. 

Different approximations are relevant to different physical situations. The projection operator 
technique has essentially the same philosophy. But with using the constraint ()4.4p in our ap- 
proach we emphasize the fundamental and central role of the Dyson equation for calculation of 
single-particle properties of many-body systems. The problem of reducing the whole hierarchy of 
equations involving higher-order GFs by a coupled nonlinear set of integro-differential equations 
connecting the single-particle GF to the self-energy operator is rather nontrivial. A characteristic 
feature of these equations is that besides the single-particle GF they involve also higher-order 
GF. The irreducible counterparts of the GFs, vertex functions, serve to identify correctly the 
self-energy as 

M = G 1 - G~ x . 

The integral form of the Dyson equation (|4.10p gives M the physical meaning of a nonlocal and 
energy-dependent effective single-particle potential. This meaning can be verified for the exact 
self-energy using the diagrammatic expansion for the causal GF. 

It is important to note that for the retarded and advanced GFs, the notion of the proper part 
M = (P) p is symbolic in natureP^EH j n a certain sense, it is possible to say that it is defined 
here by analogy with the irreducible many-particle T-matrix. Furthermore, by analogy with the 
diagrammatic technique, we can also introduce the proper part defined as a solution to the integral 
equation (|4.9p . These analogues allow us to better understand the formal structure of the Dyson 
equation for the double-time thermal GF, but only in a symbolic form . However, because of the 
identical form of the equations for GFs for all three types (advanced, retarded, and causal), we 
can convert our calculations to causal GF at each stage of calculations and, thereby, confirm the 
substantiated nature of definition (|4.9p . We therefore should speak of an analogy of the Dyson 
equation. Hereafter, we drop this stipulating, since it does not cause any misunderstanding. In a 
sense, the IGF method is a variant of the Gram-Schmidt orthogonalization procedure.^ 
It should be emphasized that the scheme presented above gives just a general idea of the IGF 
method. A more exact explanation why one should not introduce the approximation already in P, 
instead of having to work out M, is given below when working out the application of the method 
to specific problems. 

The general philosophy of the IGF method is in the separation and identification of elastic scat- 
tering effects and inelastic ones. This latter point is quite often underestimated, and both effects 
are mixed. However, as far as the right definition of quasiparticle damping is concerned, the 
separation of elastic and inelastic scattering processes is believed to be crucially important for 
many-body systems with complicated spectra and strong interaction. 

From a technical point of view, the elastic GMF renormalizations can exhibit quite a nontrivial 
structure. To obtain this structure correctly, one should construct the full GF from the complete 
algebra of relevant operators and develop a special projection procedure for higher-order GFs, in 
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accordance with a given algebra. Then a natural question arises how to select the relevant set of 
operators {Ai,^, -.-An} describing the "relevant degrees of freedom". The above consideration 
suggests an intuitive and heuristic way to the suitable procedure as arising from an infinite chain 
of equations of motion (14. 2p . Let us consider the column 

MA 

A 2 



\A n J 



where 



A 1 = A, A 2 = [A,H], A 3 = [[A,H],H},...A n = [[...[A,H}...H] 

n 

Then the most general possible Green function can be expressed as a matrix 



G = (( 



A 2 
\A n J 



a\ a\ 



At))) 



This generalized Green function describes the one-, two-, and n-particle dynamics. The equa- 
tion of motion for it includes, as a particular case, the Dyson equation for single-particle Green 
function, and the Bethe-Salpeter equation which is the equation of motion for the two-particle 
Green function and which is an analogue of the Dyson equation, etc . The corresponding reduced 
equations should be extracted from the equation of motion for the generalized Green function 
with the aid of special techniques such as the projection method and similar techniques. This 
must be a final goal towards a real understanding of the true many-body dynamics. At this 
point, it is worthwhile to underline that the above discussion is a heuristic scheme only, but not 
a straightforward recipe. The specific method of introducing the IGFs depends on the form of 
operators A n , the type of the Hamiltonian, and conditions of the problem. Here a sketchy form 
of the IGF method was presented. The aim was to introduce the general scheme and to lay the 
groundwork for generalizations. We demonstrated irp2ll2SI that the IGF method is a powerful tool 
for describing the quasiparticle excitation spectra, allowing a deeper understanding of elastic and 
inelastic quasiparticle scattering effects and the corresponding aspects of damping and finite life- 
times. In a certain sense, it provides a clear link between the equation-of-motion approach and the 
diagrammatic methods due to derivation of the Dyson equation. Moreover, due to the fact that it 
allows the approximate treatment of the self-energy effects on a final stage, it yields a systematic 
way of the construction of approximate solutions. 



5 Effective and Generalized Mean Fields 

The most common technique for studying the subject of interacting many-particle systems is 
to use the mean field theory. This approximation was especially popular in the theory of mag- 
netismP^^S To calculate the susceptibility and other characteristic functions of a system of 
localized magnetic moments, with a given interaction Hamiltonian, the approximation, termed 
the "molecular field approximation" was used widely. However, it is not an easy task to give the 
formal unified definition what the mean field is. In a sense, the mean field is the umbrella term 
for a variety of theoretical methods of reducing the many-particle problem to the single-particle 
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one. Mean field theory, that approximates the behaviour of a system by ignoring the effect of fluc- 
tuations and those spin correlations which dominate the collective properties of the ferromagnet 
usually provides a starting and estimating point only, for studying phase transitions. The mean 
field theories miss important features of the dynamics of a system. The main intention of the 
mean field theories, starting from the works of J. D. van der Waals and P. Weiss, is to take into 
account the cooperative behaviour of a large number of particles. It is well known that earlier 
theories of phase transitions based on the ideas of van der Waals and Weiss lead to predictions 
which are qualitatively at variance with results of measurements near the critical point. Other 
variants of simplified mean field theories such as the Hartree-Fock theory for electrons in atoms, 
lead to discrepancies of various kinds too. It is therefore natural to analyze the reasons for such 
drawbacks of earlier variants of the mean field theories. 

A number of effective field theories which are improved versions of the "molecular field approx- 
imation" were proposed. In our papers^K^^^ we stressed a specificity of strongly correlated 
many-particle systems on a lattice contrary to continuum (uniform) systems. The earlier concepts 
of molecular field were described in terms of a functional of mean magnetic moments (in magnetic 
terminology) or mean particle densities. The corresponding mean-field functional F[(n), {S z )] de- 
scribes the uniform mean field. Actually, the Weiss model was not based on discrete "spins" as 
is well known, but the uniformity of the mean internal field was the most essential feature of the 
model. In the modern language, one should assume that the interaction between atomic spins 
Si and its neighbors is equivalent to a mean (or molecular) field, M{ = Xo[h\ e ^ + hf" 1 ] and 
that the molecular field h\ m ^ is of the form h^ m ^ = J2i J{Rji)(Si) (above T c ). Here h ext is an 
applied conjugate field, xo is the response function, and J(Rji) is an interaction. In other words, 
the mean field approximation reduces the many-particle problem to a single-site problem in which 
a magnetic moment at any site can be either parallel or antiparallel to the total magnetic field 
composed of the applied field and the molecular field. The average interaction of % neighbors was 
taken into account only, and the fluctuations were neglected. One particular example, where the 
mean field theory works relatively well is the homogeneous structural phase transitions; in this case 
the fluctuations are confined in phase space. The next important step was made by L. Neel. He 
conjectured that the Weiss internal field might be either positive or negative in sign. In the latter 
case, he showed that below a critical temperature (Neel temperature) an ordered arrangement of 
equal numbers of oppositely directed atomic moments could be energetically favorable. This new 
magnetic structure was termed antiferromagnetism. It was conjectured that the two-sublattice 
Neel (classical) ground state is formed by local staggered internal mean fields. 
There is a number of the "correlated effective field" theories, that tend to repair the limitations 
of simplified mean field theories. The remarkable and ingenious one is the Onsager "reaction field 
approximation". He suggested that the part of the molecular field on a given dipole moment 
which comes from the reaction of neighboring molecules to the instantaneous orientation of the 
moment should not be included into the effective orienting field. This "reaction field" simply fol- 
lows the motion of the moment and thus does not favor one orientation over another (for details 
see RefePEHI). 

It is knownP^l that mean-field approximations, for example the molecular field approximation for 
a spin system, the Hartree-Fock approximation and the BCS-Bogoliubov approximation for an 
electron system are universally formulated by the Bogoliubov inequality: 

- p-HnfTret-W) < 

-p-MTveM* >) + Trei .; Hmf) l - (5-1) 

Here F is the free energy, and H m f is a "trial" or a "mean field" approximating Hamiltonian. 
This inequality gives the upper bound of the free energy of a many-body system. It is important 
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to emphasize that the BCS-Bogoliubov theory of superconductivity^"^^ was formulated on 
the basis of a trial Hamiltonian which consists of a quadratic form of creation and annihilation 
operators, including "anomalous" (off-diagonal) averages. The functional of the mean field (for 
the superconducting single-band Hubbard model) is of the following formP3 



1 M-A.) -<4^> 1 



The "anomalous" off-diagonal terms fix the relevant BCS-Bogoliubov vacuum and select the ap- 
propriate set of solutions. From the point of view of quantum many-body theory, the problem of 
adequate introduction of mean fields for system of many interacting particles can be most consis- 
tently investigated in the framework of the IGF method. A correct calculation of the quasiparticle 
spectra and their damping, particularly, for systems with a complicated spectrum and strong in- 
teraction^ reveals, that the generalized mean fields can have very complicated structure which 
cannot be described by a functional of the mean-particle density. 

To illustrate the actual distinction of description of the generalized mean field in the equation-of- 
motion method for the double-time Green functions, let us compare the two approaches, namely, 
that of Tyablikov-^ and of C alien P*^ We shall consider the Green function ((S + \S~)} for the 
isotropic Heisenberg model 

H = ~J2j(i-j)S i S j . (5.3) 

The equation of motion for the spin Green function is of the form 

u((S+\Sr)) u = (5.4) 
2{S z )5 l3 + Y J J(i-g){{StS z g - Spt\Sj)) w . 

9 

The Tyablikov decoupling expresses the second-order Green function in terms of the first (initial) 
Green function: 

((S+S z g \Sr)) = (S z )((S+\Sr)). (5.5) 
This approximation is an RPA-type; it does not lead to the damping of spin wave excitations 

E(q) = Y,J( i -9)(S z )ex P [i(R i -R g )q\=2(S z )(Jo-J q ). (5.6) 

9 

The reason for this is rather transparent. This decoupling does not take into account the inelastic 
magnon-magnon scattering processes. In a sense, the Tyablikov approximation consists of approx- 
imating the commutation relations of spin operators to the extent of replacing the commutation 
relation = 2S?5 tJ by [S+,Sr}_ = 2(S z )5 ij . 

CallerPSl has proposed an improved decoupling approximation in the method of Tyablikov in the 
following form: 

(S z )((SpB)) - a(S g Sp((S+\B)). 

Here < a < 1. To clarify this point, it should be reminded that for spin 1/2 ( the procedure was 
generalized by Callen to an arbitrary spin), the spin operator S z can be written as Sg = S — S~S+ 
or Sg = ^(S^S~ — 5 ~ 5+ ) . It is easy to show that 

Sg =aS H — S g S g — Sg ^ . 
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The operator S g S£ represents the deviation of (S z ) from S. In the low-temperature region, this 
deviation is small, and a ~ 1. Similarly, the operator ^(S^S~ — S~S£) represents the deviation 
of (S z ) from 0. Thus, when (S z ) approaches to zero, one can expect that a ~ 0. Thus, in this 
way, it is possible to obtain a correction to the Tyablikov decoupling with either a positive or 
negative sign, or no correction at all, or any intermediate value, depending on the choice of a. 
The above Callen arguments are not rigorous , for, although the difference in the operators S + S~ 
and S~S + is small if (S z ) ~ 0, each operator makes a contribution of the order of S, and it is 
each operator which is treated approximately, not the difference. There are some other drawbacks 
of the Callen decoupling scheme. Nevertheless, the Callen decoupling was the first conceptual 
attempt to introduce the interpolation decoupling procedure. Let us note that the choice of a = 
over the entire temperature range is just the Tyablikov decoupling (|5.5p . 
The energy spectrum for the Callen decoupling is given by 

E(q) = 2(S Z )((J - J q ) + ^£[J(fc) - J{k - q)]N{E(k))). (5.7) 

k 

Here N{E{k)) is the Bose distribution function N(E(k)) = [exp(E(k)f3) -l] -1 . This is an implicit 
equation for N(E(k)), involving the unknown quantity {S z ) . For the latter an additional equation 
is givenPSl Thus, both these equations constitute a set of coupled equations which must be solved 
self-consistently for (S z ). This formulation of the Callen decoupling scheme displays explicitly 
the tendency of the improved description of the mean field. In a sense, it is possible to say 
that the Callen work dates really the idea of the generalized mean field within the equation-of- 
motion method for double-time GFs, however, in a semi-intuitive form. The next essential steps 
were made by PlakidcP^ for the Heisenberg ferromagnet and by Kuzemsky^ for the Hubbard 
model. Later many approximate schemes for decoupling the hierarchy of equations for GF were 
proposed, improving the Tyablikov's and Callen's decouplings. Various approaches generalizing 
the random phase's approximation in the ferromagnetism theory for wide ranges of temperature 
were considered in paper by Czachor and HolasP^ As was mentioned above, the correct definition 
of generalized mean fields depends on the condition of the problem, the strength of interaction, the 
choice of relevant operators, and on the symmetry requirements. The most important conclusion 
to be drawn from the present consideration is that the GMF, in principle, can have quite a 
nontrivial structure and cannot be reduced to the mean-density functional only. 

6 Quasiaverages and Irreducible Green Functions Method 

In condensed matter physics, the symmetry is important in classifying different phases and under- 
standing the phase transitions between them. There is an important distinction between the case 
where the broken symmetry is continuous (e.g. translation, rotation, gauge invariance) or discrete 
(e.g. inversion, time reversal symmetry) P2 The Goldstone theorem states that when a continuous 
symmetry is spontaneously broken and the interactions are short ranged a collective mode (excita- 
tion) exists with a gapless energy spectrum (i.e. the energy dispersion curve starts at zero energy 
and is continuous). Acoustical phonons in a crystal are prime examples of such so-called gapless 
Goldstone modes. Other examples are the Bogoliubov sound modes in (charge neutral) Bose 
condensates32HSl a nd spin waves (magnons) in ferro- and antiferromagnets. N. N. Bogoliubov and 
then Y. Nambu in in their works shown that the general features of superconductivity are in fact 
model independent consequences of the spontaneous breakdown of electromagnetic gauge invari- 
ance. It is important to emphasize that the BCS-Bogoliubov theory of superconductivitj^^IES] 
was formulated on the basis of a trial Hamiltonian which consists of a quadratic form of creation 
and annihilation operators, including "anomalous" (off-diagonal) averages.^ 3 The strong-coupling 
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BCS-Bogoliubov theory of superconductivity was formulated for the Hubbard model in the local- 
ized Wannier representation in RefsF^®"^ Therefore, instead of the algebra of the normal state's 
operator a^, a\ a and rii a , for description of superconducting states, one has to use a more general 
algebra, which includes the operators di a , a\ a , rii a and ai a ai- a , a\ a a\_ a . The relevant generalized 
one-electron Green function will have the following formPSl^^Ql 

(Gu G 12 \ ( {{aia\a) a )) {{a ia \aj- a ))\ 

As it was discussed in Refs.pM^the off-diagonal (anomalous) entries of the above matrix select 
the vacuum state of the system in the BCS-Bogoliubov form, and they are responsible for the 
presence of anomalous averages. For treating the problem we follow the general scheme of the 
irreducible Green functions methodJ22E^ In this approach we start from the equation of motion 
for the Green function Gij(ui) (normal and anomalous components) 

y^juSjj - tij)((a ja \al a )) = 5u, + (6.2) 
j 

U {{a i(T ni- a \al /a . crUn\ala))i 

nj 

X>^ + %)((a]->L>>= ( 6 - 3 ) 
i 

- U {{ a \-a n i°\ a \>o)) + ^2 V M(a]_ a u n \a], a }}. 

nj 

The irreducible Green functions are introduced by definition 

(^((oiaot^oi-alot^)),,,) = {{a ia a\_ a ai^ a \a\, a )) ^ - (6.4) 
-(ni- a )Gn + (aj (T a i _ (7 )((aJ_ (T |a t T , (T )) w , 

( M ((4* a i«4~*\4*))u) = ((4x a ^4->L>>" - 
-{n ia )G 2 i + (ata|_ (J )((a i(T |a l t , (T )) tJ . 

The self-consistent system of superconductivity equations follows from the Dyson equatiorP^^EB 

G tl ,(u) = G% (to) +Y J G Q ij {u)M jf {u)G j , i ,{u>). (6.5) 
jj' 

The mass operator Mjj>(uj) describes the processes of inelastic electron scattering on lattice vibra- 
tions. The elastic processes are described by the quantity (|5.2|) . Thus the "anomalous" off-diagonal 
terms fix the relevant BCS-Bogoliubov vacuum and select the appropriate set of solutions. The 
functional of the generalized mean field for the superconducting single-band Hubbard model is of 
the form Yf a . A remark about the BCS-Bogoliubov mean-field approach is instructive. Speaking 
in physical terms, this theory involves a condensation correctly, in spite that such a condensation 
cannot be obtained by an expansion in the effective interaction between electrons. Other mean 
field theories, e.g. the Weiss molecular field theory and the van der Waals theory of the liquid-gas 
transition are much less reliable. The reason why a mean-field theory of the superconductivity 
in the BCS-Bogoliubov form is successful would appear to be that the main correlations in metal 
are governed by the extreme degeneracy of the electron gas. The correlations due to the pair 
condensation, although they have dramatic effects, are weak (at least in the ordinary supercon- 
ductors) in comparison with the typical electron energies, and may be treated in an average way 
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with a reasonable accuracy. It should be emphasized that the high-temperature superconductors, 
discovered two decades ago, motivated an intensification of research in superconductivity, not only 
because applications are promising, but because they also represent a new state of matter that 
breaks certain fundamental symmetries. These are the broken symmetries of gauge (superconduc- 
tivity), reflection (d-wave superconducting order parameter), and time-reversal (ferromagnetism) . 
Superconductivity and antiferromagnetism are both the spontaneously broken symmetries. The 
question of symmetry breaking within the localized and band models of antiferromagnets was 
studied by the author of this work in Refs! 22 ^ 32 ^ 3 It has been found there that the concept of 
spontaneous symmetry breaking in the band model of magnetism is much more complicated than 
in the localized model. In the framework of the band model of magnetism one has to additionally 
consider the so called anomalous propagators of the form 

FM:G /m ~((a to ;4j), 
AFM : G afm ~ ((a k+Q(T ;al +Q , a ,)). 

In the case of the band antiferromagnet the ground state of the system corresponds to a spin- 
density wave (SDW), where a particle scattered on the internal inhomogeneous periodic field gains 
the momentum Q — Q' and changes its spin: a — >• & ' . The long-range order parameters are defined 
as follows: 

FM : m = l/N ^ a L a h-a), (6.6) 

ka 

AFM : M Q = Y,(4* a k+Q-*)- (6.7) 

ka 

It is important to stress, that the long-range order parameters here are functionals of the internal 
field, which in turn is a function of the order parameter. Thus, in the cases of rotation and trans- 
lation invariant Hamiltonians of band ferro- and antiferromagnetics one has to add the following 
infinitesimal sources removing the degeneracy: 

FM:u^ B H x J2<it a k-a, (6.8) 

ka 

AFM : ufi B H 4a a k+Q-a- (6.9) 
kQ 

Here, v — )■ after the usual in statistical mechanics infinite-volume limit V — > oo. The ground 
state in the form of a spin-density wave was obtained for the first time by Overhauser. There, 
the vector Q is a measure of inhomogeneity or translation symmetry breaking in the system. The 
analysis performed by various authors showed that the antiferromagnetic and more complicated 
states (for instance, ferrimagnetic) can be described in the framework of a generalized mean-field 
approximation.^ In doing that we have to take into account both the normal averages {a\ a ai a ) 
and the anomalous averages (al a ai^ a }. It is clear that the anomalous terms break the original 
rotational symmetry of the Hubbard Hamiltonian. Thus, the generalized mean-field's approxi- 
mation for the antiferromagnet has the following formp3 m- a ai a ~ (nj_ (7 )aj - — (a|_ cr aj -)aj_ cr . A 
self-consistent theory of band antiferromagnetisrrPH was developed by the author of this work us- 
ing the method of the irreducible Green functions ! 22 H 38 I The following definition of the irreducible 
Green functions was used: 

{{ a k+pa a \>-\-q-a a q-a\ a \ a )) 'ui = ( (o-k+pa a p+q-a a q~a \ a ka) ) w ~~ 

Sp,o(n q -a)G ka - (ak+pa<4,+q-a)((aq~aW ka ))uj- {6.10) 
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The algebra of relevant operators must be chosen as follows ((cw, aj CT , n^, a\ a ai- a ). The corre- 
sponding initial GF will have the following matrix structure 



Qafm 



\a) a )) ((a la \a]_ a )) 



J) {{«-*\a}_J) 



The off-diagonal terms select the vacuum state of the band's antiferromagnet in the form of a spin- 
density wave. With this definition, one introduces the so-called anomalous (off-diagonal) Green 
functions which fix the relevant vacuum and select the proper symmetry broken solutions. The 
theory of the itinerant antiferromagnetisirpH was formulated by using sophisticated arguments of 
the irreducible Green functions method in complete analogy with our description of the Heisenberg 
antiferromagnet at finite temperatures.^ For the two-sublattice antiferromagnet we used the 
matrix Green function of the form 

G(k ^-{((s+\s Zka )) ((s+\s Zkb )))- t 6 - 11 ) 

Here, the Green functions on the main diagonal are the usual or normal Green functions, while the 
off-diagonal Green functions describe contributions from the so-called anomalous terms, analogous 
to the anomalous terms in the BCS-Bogoliubov superconductivity theory. The anomalous (or off- 
diagonal) average values in this case select the vacuum state of the system precisely in the form 
of the two-sublattice Neel stated 



7 Conclusions 

In the present work we shown that the development and improvement of the methods of quantum 
statistical mechanics still remains quite an important direction of research. In particular, the 
Bogoliubov's method of quasiaverages gives the deep foundation and clarification of the concept of 
broken symmetry.^ It makes the emphasis on the notion of a degeneracy and plays an important 
role in equilibrium statistical mechanics of many-particle systems. According to that concept, 
infinitely small perturbations can trigger macroscopic responses in the system if they break some 
symmetry and remove the related degeneracy (or quasidegeneracy) of the equilibrium state. As 
a result, they can produce macroscopic effects even when the perturbation magnitude is tend to 
zero, provided that happens after passing to the thermodynamic limit. 

We have discussed the theory of the correlation effects for many-particle interacting systems 
using the ideas of quasiaverages for interacting electron and spin systems on a lattice. The 
workable and self-consistent irreducible Green functions approach to the decoupling problem for 
the equation-of-motion method for double-time temperature Green functions has been presented. 
The main advantage of the formalism consists in the clear separation of the elastic scattering 
corrections (generalized mean fields) and inelastic scattering effects (damping and finite lifetimes). 
These effects could be self-consistently incorporated in a general and compact manner. Using 
the IGF method, it is possible to obtain a closed self-consistent set of equations determining the 
relevant Green functions and self-energy. These equations give a general microscopic description of 
correlation effects. Moreover, this approach gives the workable scheme for the definition of relevant 
generalized mean fields written in terms of appropriate correlators. This picture of interacting 
many-particle systems on a lattice is far richer and gives more possibilities for the analysis of 
phenomena which can actually take place. In this sense the approach we described produces more 
advanced physical picture of the quasiparticle many-body dynamics. 
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Abstract 

The development and applications of the method of quasiaverages developed by N. N. 
Bogoliubov to quantum statistical physics and to quantum solid state theory and, in particular, 
to quantum theory of magnetism, were analyzed. The problem of finding the ferromagnetic, 
antiferromagnetic and superconducting symmetry broken solutions of the correlated lattice 
fermion models was discussed within the irreducible Green functions method. A unified scheme 
for the construction of generalized mean fields (elastic scattering corrections) and self-energy 
(inelastic scattering) in terms of the Dyson equation was generalized in order to include the 
source fields. The interrelation of the Bogoliubov's idea of quasiaverages and the concepts of 
symmetry breaking and quantum protectorate was discussed briefly in the context of quantum 
statistical physics. The idea of quantum protectorate reveals the essential difference in the 
behavior of the complex many-body systems at the low-energy and high-energy scales. It was 
shown that the role of symmetry (and the breaking of symmetries) in combination with the 
degeneracy of the system was reanalyzed and essentially clarified within the framework of the 
method of quasiaverages. The complementary notion of quantum protectorate might provide 
distinctive signatures and good criteria for a hierarchy of energy scales and the appropriate 
emergent behavior. 
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1 Introduction 



It is well known that symmetry principles play a crucial role in physicsP^ The theory of symmetry 
is a basic tool for understanding and formulating the fundamental notions of physics P According 
to F. WilczekS "the primary goal of fundamental physics is to discover profound concepts that 
illuminate our understanding of nature" . It is known that symmetry is a driving force in the shap- 
ing of physical theory; moreover, the primary lesson of physics of last century is that the secret 
of nature is symmetry. Every symmetry leads to a conservation law; the well known examples 
are the conservation of energy, momentum and electrical charge. A variety of other conservation 
laws can be deduced from symmetry or invariance properties of the corresponding Lagrangian or 
Hamiltonian of the system. According to Noether theorem, every continuous symmetry transfor- 
mation under which the Lagrangian of a given system remains invariant implies the existence of 
a conserved function. 2 Many fundamental laws of physics in addition to their detailed features 
possess various symmetry properties. These symmetry properties lead to certain constraints and 
regularities on the possible properties of matter. Thus the principles of symmetries belong to the 
underlying principles of physics. 

It is known that when the Hamiltonian of a system is invariant under a symmetry operation, but 
the ground state is not, the symmetry of the system can be spontaneously brokenP Symmetry 
breaking is termed spontaneous when there is no explicit term in a Lagrangian which manifestly 
breaks the symmetry. Symmetries and breaking of symmetries play an important role in statistical 
physics, quantum field theory, physics of elementary particles, etcP^ 

In physics, spontaneous symmetry breaking occurs when a system that is symmetric with respect 
to some symmetry group goes into a vacuum state that is not symmetric. When that happens, the 
system no longer appears to behave in a symmetric manner. It is a phenomenon that naturally 
occurs in many situations. The mechanism of spontaneous symmetry breaking is usually under- 
stood as the mechanism responsible for the occurrence of asymmetric states in quantum systems 
in the thermodynamic limit and is used in various field of quantum physicsP The intriguing 
mechanism of spontaneous symmetry breaking is a unifying concept that lie at the basis of most 
of the recent developments in theoretical physics, from statistical mechanics to many-body theory 
and to elementary particles theoryP^ 

It should be stressed that symmetry implies degeneracy. The greater the symmetry, the greater 
the degeneracy. The study of the degeneracy of the energy levels plays a very important role in 
quantum physics. It is of importance to emphasize that when spontaneous symmetry breaking 
takes place, the ground state of the system is degenerate. Substantial progress in the under- 
standing of the broken symmetry concept was connected with Bogoliubov's fundamental ideas on 
quasiaverages Studies of degenerate systems led Bogoliubov in 1960-61 to the formulation 
of the method of quasiaverages. This method has proved to be a universal tool for systems 
whose ground states become unstable under small perturbations. Thus the role of symmetry (and 
the breaking of symmetries) in combination with the degeneracy of the system was reanalyzed 
and essentially clarified by N. N. Bogoliubov in 1960-1961. He invented and formulated a powerful 
innovative idea of quasiaverages in statistical mechanics JSEMffllMISllIS The very elegant work of 
N. N. Bogoliubov^ has been of great importance for a deeper understanding of phase transitions, 
superfluidity and superconductivity, quantum theory of magnetism^ and other fields of equilib- 
rium and nonequilibrium statistical rnechanics.^^HIllElIIElliQHlS] ^he concept of quasiaverages is 
indirectly related to the theory of phase transition. The instability of thermodynamic averages 
with respect to perturbations of the Hamiltonian by a breaking of the invariance with respect 
to a certain group of transformations means that in the system transition to an extremal state 
occurs. The mathematical apparatus of the method of quasiaverages includes the Bogoliubov 
theorerr F^^ ^^ on singularities of type 1/g 2 and the Bogoliubov inequality for Green and cor- 
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relation functions as a direct consequence of the method. It includes algorithms for establishing 
non-trivial estimates for equilibrium quasiaverages, enabling one to study the problem of ordering 
in statistical systems and to elucidate the structure of the energy spectrum of the underlying 
excited states. Thus the Bogoliubov's idea of quasiaverages is an essential conceptual advance of 
modern physics. 

Many-particle systems where the interaction is strong have often complicated behavior, and re- 
quire nonperturbative approaches to treat their properties. Such situations are often arise in 
condensed matter systems. Electrical, magnetic and mechanical properties of materials are emer- 
gent collective behaviors of the underlying quantum mechanics of their electrons and constituent 
atoms. A principal aim of solid state physics and materials science is to elucidate this emergence. 
A full achievement of this goal would imply the ability to engineer a material that is optimum 
for any particular application. The current understanding of electrons in solids uses simplified 
but workable picture known as the Fermi liquid theory. This theory explains why electrons in 
solids can often be described in a simplified manner which appears to ignore the large repulsive 
forces that electrons are known to exert on one another. There is a growing appreciation that this 
theory probably fails for entire classes of possibly useful materials and there is the suspicion that 
the failure has to do with unresolved competition between different possible emergent behaviors. 
It is appropriate to note here that the emergent properties of matter were analyzed and discussed 
by R. Laughlin and D. Pines^EID f r0 m a general point of view (see also RefP^J). They introduced 
a unifying idea of quantum protectorate. This concept belongs also to the underlying principles 
of physics. The idea of quantum protectorate reveals the essential difference in the behavior of 
the complex many-body systems at the low-energy and high-energy scales. The existence of two 
scales, low-energy and high-energy, in the description of physical phenomena is used in physics, 
explicitly or implicitly. 

From the other hand, it was recognized for many years that the strong correlation in solids exist 
between the motions of various particles (electrons and ions, i.e. the fermion and boson degrees 
of freedom) which arise from the Coulomb forces. The most interesting objects are metals and 
their compounds. They are invariant under the translation group of a crystal lattice and have 
lattice vibrations as well as electron degrees of freedom. There are many evidences for the im- 
portance of many-body effects in these systems. Within the semi-phenomenological theory it was 
suggested that the low-lying excited states of an interacting Fermi gas can be described in terms 
of a set of "independent quasiparticles". However, this was a phenomenological approach and 
did not reveal the nature of relevant interactions. An alternative way of viewing quasiparticles, 
more general and consistent, is through the Green function scheme of many-body theory IISIEZIEBIES] 
It becomes clear that only a thorough experimental and theoretical investigation of quasiparti- 
cle many-body dynamics of the many particle systems can provide the answer on the relevant 
microscopic pictured In our works, we discussed the microscopic view of a dynamic behavior 
of various interacting many-body systems on a latticeP^SBS] A comprehensive description of 
transition and rare-earth metals and alloys and other materials (as well as efficient predictions of 
properties of new materials) is possible only in those cases, when there is an adequate quantum- 
statistical theory based on the information about the electron and crystalline structures. The 
main theoretical problem of this direction of research, which is the essence of the quantum theory 
of magnetism, is investigations and improvements of quantum-statistical models describing the 
behavior of the complex compounds and materials in order to take into account the main features 
of their electronic structure, namely, their dual "band-atomic" nature.^ The construction of a 
consistent theory explaining the electronic structure of these substances encounters serious diffi- 
culties when trying to describe the collectivization-localization duality in the behavior of electrons. 
This problem appears to be extremely important, since its solution gives us a key to understanding 
magnetic, electronic, and other properties of this diverse group of substances. The author of the 
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present paper investigated the suitability of the basic models with strong electron correlations and 
with a complex spectrum for an adequate and correct description of the dual character of electron 
statespZl A universal mathematical formalism was developed for this investigationPSEH n takes 
into account the main features of the electronic structure and allows one to describe the true 
quasiparticle spectrum, as well as the appearance of the magnetically ordered, superconducting, 
and dielectric (or semiconducting) states. With a few exceptions, diverse physical phenomena 
observed in compounds and alloys of transition and rare-earth metals, cannot be explained in 
the framework of the mean- field approximation, which overestimates the role of inter-electron 
correlations in computations of their static and dynamic characteristics. The circle of questions 
without a precise and definitive answer, so far, includes such extremely important (not only from 
a theoretical, but also from a practical point of view) problems as the adequate description of 
quasiparticle dynamics for quantum-statistical models in a wide range of their parameter values. 
The source of difficulties here lies not only in the complexity of calculations of certain dynamic 
properties (such as, the density of states, electrical conductivity, susceptibility, electron-phonon 
spectral function, the inelastic scattering cross section for slow neutrons), but also in the ab- 
sence of a well-developed method for a consistent quantum-statistical analysis of a many-particle 
interaction in such systems. A self-consistent field approach was used in the papers^E^SS] f or 
description of various dynamic characteristics of strongly correlated electronic systems. It allows 
one to consistently and quite compactly compute quasiparticle spectra for many-particle systems 
with strong interaction taking into account damping effects. The correlation effects and quasipar- 
ticle damping are the determining factors in analysis of the normal properties of high-temperature 
superconductors, heavy fermion compounds, etc. We also formulated a general scheme for a the- 
oretical description of electronic properties of many-particle systems taking into account strong 
inter-electron correlations! 22 ^ 38 ^ The scheme is a synthesis of the method of two-time tempera- 
ture Green's functionaries] anc i th e diagram technique. An important feature of this approach is 
a clear-cut separation of the elastic and inelastic scattering processes in many-particle systems 
(which is a highly nontrivial task for strongly correlated systems). As a result, one can construct a 
correct basic approximation in terms of generalized mean fields (the elastic scattering corrections), 
which allows one to describe magnetically ordered or superconducting states of the system. The 
residual correlation effects, which are the source of quasiparticle damping, are described in terms 
of the Dyson equation with a formally exact representation for the mass operator. 
In the present paper we will discuss some applications of the symmetry principles to quantum 
and statistical physics and quantum solid state theory in the light of our results on quasiparticle 
many-body dynamics. 

2 Bogoliubov's Quasiaverages in Statistical Mechanics 

In the work of N. N. Bogoliubov "Quasiaverages in Problems of Statistical Mechanics" the innova- 
tive notion of quasiavereg^^ was introduced and applied to various problem of statistical physics. 
In particular, quasiaverages of Green's functions constructed from ordinary averages, degeneration 
of statistical equilibrium states, principle of weakened correlations, and particle pair states were 
considered. In this framework the l/t/ 2 -type properties in the theory of the superfluidity of Bose 
and Fermi systems, the properties of basic Green functions for a Bose system in the presence of 
condensate, and a model with separated condensate were analyzed. 

The method of quasiaverages is a constructive workable scheme for studying systems with spon- 
taneous symmetry breakdown. A quasiaverage is a thermodynamic (in statistical mechanics) or 
vacuum (in quantum field theory) average of dynamical quantities in a specially modified averag- 
ing procedure, enabling one to take into account the effects of the influence of state degeneracy of 
the system. The method gives the so-called macro-objectivation of the degeneracy in the domain 
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of quantum statistical mechanics and in quantum physics. In statistical mechanics, under sponta- 
neous symmetry breakdown one can, by using the method of quasiaverages, describe macroscopic 
observable within the framework of the microscopic approach. 

In considering problems of findings the eigenfunctions in quantum mechanics it is well known 
that the theory of perturbations should be modified substantially for the degenerate systems. In 
the problems of statistical mechanics we have always the degenerate case due to existence of the 
additive conservation laws. The traditional approach to quantum statistical mechanics^ is based 
on the unique canonical quantization of classical Hamiltonians for systems with finitely many 
degrees of freedom together with the ensemble averaging in terms of traces involving a statistical 
operator p. For an operator A corresponding to some physical quantity A the average value of A 
will be given as 

(A) H = TrpA; p = exp"^ /Trexp~ PH , (2-1) 

where H is the Hamiltonian of the system, j3 = 1/kT is the reciprocal of the temperature. 
The core of the problem lies in establishing the existence of a thermodynamic limit (such as 
N/V = const, V — > oo, N = number of degrees of freedom, V = volume) and its evaluation 
for the quantities of interest. Thus in the statistical mechanics the average (A) of any dynamical 
quantity A is defined in a single- valued way. In the situations with degeneracy the specific problems 
appear. In quantum mechanics, if two linearly independent state vectors (wavefunctions in the 
Schroedinger picture) have the same energy, there is a degeneracy. In this case more than one 
independent state of the system corresponds to a single energy level. If the statistical equilibrium 
state of the system possesses lower symmetry than the Hamiltonian of the system (i.e. the situation 
with the symmetry breakdown), then it is necessary to supplement the averaging procedure ()2. 1|) 
by a rule forbidding irrelevant averaging over the values of macroscopic quantities considered 
for which a change is not accompanied by a change in energy. This is achieved by introducing 
quasiaverages, that is, averages over the Hamiltonian H v g supplemented by infinitesimally-small 
terms that violate the additive conservations laws H v g = H + u{e ■ M), (y — >■ 0). Thermodynamic 
averaging may turn out to be unstable with respect to such a change of the original Hamiltonian, 
which is another indication of degeneracy of the equilibrium state. According to Bogoliubovp^ 
the quasiaverage of a dynamical quantity A for the system with the Hamiltonian H u g is defined 
as the limit 

< A^=lim(A) u g, (2.2) 

v — ^0 

where {A) u g denotes the ordinary average taken over the Hamiltonian H u g, containing the small 
symmetry-breaking terms introduced by the inclusion parameter v, which vanish as v — > after 
passage to the thermodynamic limit V — > oo. It is important to note that in this equation limits 
cannot be interchanged. Thus the existence of degeneracy is reflected directly in the quasiaverages 
by their dependence upon the arbitrary unit vector e. It is also clear that 

(A) = J <A> de. (2.3) 

According to definition (12.3p . the ordinary thermodynamic average is obtained by extra averaging 
of the quasiaverage over the symmetry-breaking group. Thus to describe the case of a degenerate 
state of statistical equilibrium quasiaverages are more convenient, more physical, than ordinary 
averages.^ The latter are the same quasiaverages only averaged over all the directions e. 
It is necessary to stress, that the starting point for Bogoliubov's workP^ was an investigation of 
additive conservation laws and selection rules, continuing and developing the approach by P. Curie 
for derivation of selection rules for physical effects (see also RefP). Bogoliubov demonstrated that 
in the cases when the state of statistical equilibrium is degenerate, as in the case of the Heisenberg 
ferromagnet, one can remove the degeneracy of equilibrium states with respect to the group of 
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spin rotations by including in the Hamiltonian H an additional noninvariant term vM z V with an 
infinitely small v. For the Heisenberg ferromagnet the ordinary averages must be invariant with 
regard to the spin rotation group. The corresponding quasiaverages possess only the property of 
covariance. Thus the quasiaverages do not follow the same selection rules as those which govern 
ordinary averages, due to their invariance with regard to the spin rotation group. It is clear that 
that the unit vector e, i.e., the direction of the magnetization M vector, characterizes the degen- 
eracy of the considered state of statistical equilibrium. In order to remove the degeneracy one 
should fix the direction of the unit vector e. It can be chosen to be along the z direction. Then 
all the quasiaverages will be the definite numbers. This is the kind that one usually deals with in 
the theory of ferromagnetism. 

The value of the quasi-average (12. 2D may depend on the concrete structure of the additional term 
AH = H u — H, if the dynamical quantity to be averaged is not invariant with respect to the 
symmetry group of the original Hamiltonian H. For a degenerate state the limit of ordinary 
averages (|2.3p as the inclusion parameters v of the sources tend to zero in an arbitrary fashion, 
may not exist. For a complete definition of quasiaverages it is necessary to indicate the manner 
in which these parameters tend to zero in order to ensure convergence.^ On the other hand, in 
order to remove degeneracy it suffices, in the construction of H, to violate only those additive 
conservation laws whose switching lead to instability of the ordinary average. Thus in terms of 
quasiaverages the selection rules for the correlation functions^^ that are not relevant are those 
that are restricted by these conservation laws. 

By using H u , we define the state uj{A) = (A) u and then let v tend to zero (after passing to the 
thermodynamic limit) PEUE] if a ll 

averages uj (A) get infinitely small increments under infinitely 
small perturbations v, this means that the state of statistical equilibrium under consideration is 
nondegeneratePHE' However, if some states have finite increments as v — > 0, then the state is 
degenerate. In this case, instead of ordinary averages {A)h, one should introduce the quasiaver- 
ages (I2.2h . for which the usual selection rules do not hold. 

The method of quasiaverages is directly related to the principle weakening of the correlatiorfl^^^ 
in many-particle systems. According to this principle, the notion of the weakening of the corre- 
lation, known in statistical mechanics pSHHEHl i n the- case of state degeneracy must be interpreted 
in the sense of the quasiaverages P 

The quasiaverages may be obtained from the ordinary averages by using the cluster property 
which was formulated by Bogoliubov.H' This was first done when deriving the Boltzmann equa- 
tions from the chain of equations for distribution functions, and in the investigation of the model 
Hamiltonian in the theory of superconductivityP^^^^ To demonstrate this let us consider 
averages (quasiaverages) of the form 

F(ti,x 1) ...t n ,x n ) = {...&(t 1 ,xi)...-$(tj,x j )...), (2.4) 

where the number of creation operators ^ may be not equal to the number of annihilation opera- 
tors We fix times and split the arguments (t\, x\, . . . t n ,x n ) into several clusters (. . . , t a ,x a , ...),..., 
(. . . ,tp,xp, . . .). Then it is reasonably to assume that the distances between all clusters \x a — xp\ 
tend to infinity. Then, according to the cluster property, the average value (j2.4|) tends to the prod- 
uct of averages of collections of operators with the arguments (. . . , t a ,x a , ...),..., (. . . , tp,xp, . . .) 

lim F(ti,xi,...t n ,x n ) = F(...,t a ,x a ,...)...F(...,t/3,Xi3,...). (2.5) 

Xg |— too 

For equilibrium states with small densities and short-range potential, the validity of this property 
can be proved.^ For the general case, the validity of the cluster property has not yet been proved. 
Bogoliubov formulated it not only for ordinary averages but also for quasiaverages, i.e., for anoma- 
lous averages, too. It works for many important models, including the models of superfluidity^^ 
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and superconductivitjPESEEl ( see also RefsPHSl). 

To illustrate this statement consider Bogoliubov's theory of a Bose-system with separated con- 
densate, which is given by the Hamiltoniarl^"^ 

H A = [ yU x )(- — )V(x)dx - ii ! &(x)V(x)dx (2.6) 
Ja 2m Ja 

+7; I ^{ x i)^H x 2)<£>(xi ~ x 2 )^{x 2 )^{xi)dxidx 2 . 

2 JA2 

This Hamiltonian can be written also in the following form 

H A = H + H 1= I y\ q )(-A)y( q )dq (2.7) 



+ - / y\q)&(q')<5>(q-q')^(q')y(q)dqdq'. 
' I A 2 
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Here, ^(q), and W(q) are the operators of annihilation and creation of bosons. They satisfy the 
canonical commutation relations 

[*(<?), = Kq - <?'); *(</)] = *V)] = 0. (2.8) 

The system of bosons is contained in the cube A with the edge L and volume V. It was assumed 
that it satisfies periodic boundary conditions and the potential $((?) is spherically symmetric and 
proportional to the small parameter. It was also assumed that, at temperature zero, a certain 
macroscopic number of particles having a nonzero density is situated in the state with momentum 
zero. The operators ^(q), and ^(q) can be represented in the form 

= a /W; &(q) = 4/W, (2.9) 

where ciq and are the operators of annihilation and creation of particles with momentum zero. 
To explain the phenomenon of superfluidity, one should calculate the spectrum of the Hamilto- 
nian, which is quite a difficult problem. Bogoliubov suggested the idea of approximate calculation 
of the spectrum of the ground state and its elementary excitations based on the physical nature 
of superfluidity. His idea consists of a few assumptions. The main assumption is that at temper- 
ature zero the macroscopic number of particles (with nonzero density) has the momentum zero. 
Therefore, in the thermodynamic limit, the operators ao/W and a^/VV commute 

lim 

y-»oo 



a /W,al/W]=^^0 (2.10) 



and are c-numbers. Hence, the operator of the number of particles No = a^ao is a c-number, 
too. It is worth noting that the Hamiltonian (|2.7h is invariant under the gauge transformation 
dk = exp(iip)dk, at = exp(— i<p)a\, where cp is an arbitrary real number. Therefore, the averages 
(ao/y/V) and (o,q/\/V) must vanish. But this contradicts to the assumption that clq/VV and 
a /\/V must become c-numbers in the thermodynamic limit. In addition it must be taken into 
account that a^ao/V = Nq/V ^ which implies th&tao/VV = No ex.p(ia)/W ^ and a\/y/V = 
No exp(— ia)/\/V ^ 0, where a is an arbitrary real number. This contradiction may be overcome 
if we assume that the eigenstates of the Hamiltonian are degenerate and not invariant under gauge 
transformations, i.e., that a spontaneous breaking of symmetry takes place. 

Thus the averages {ao/y/V) and (aJ/VF), which are nonzero under spontaneously broken gauge 
invariance, are called anomalous averages or quasiaverages. This innovative idea of Bogoliubov 
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penetrate deeply into the modern quantum physics. The systems with spontaneously broken 
symmetry are studied by use of the transformation of the operators of the form 

= a /W + 8(qy, &(q) = al/W + 0*(q), (2.11) 

where ao/VV and a\/\/V are the numbers first introduced by Bogoliubov in 1947 in his inves- 
tigation of the phenomenon of superfluidityP^^^ The main conclusion was made that for the 
systems with spontaneously broken symmetry, the quasiaverages should be studied instead of 
the ordinary averages. It turns out that the long-range order appears not only in the system of 
Bose-particles but also in all systems with spontaneously broken symmetry. Bogoliubov's papers 
outlined above anticipated the methods of investigation of systems with spontaneously broken 
symmetry for many years. 

As mentioned above, in order to explain the phenomenon of superfluidity, Bogoliubov assumed 
that the operators ao/VV and a^/VV become c-numbers in the thermodynamic limit. This state- 
ment was rigorously proved in the papers by Bogolyubov and some other authors. Bogolyubov's 
proof was based on the study of the equations for two-time Green's functions^ and on the as- 
sumption that the cluster property holds. It was proved that the solutions of equations for Green's 
functions for the system with Hamiltonian (|2.7p coincide with the solutions of the equations for 
the system with the same Hamiltonian in which the operators ao/VV and a /W are replaced by 
numbers. These numbers should be determined from the condition of minimum for free energy. 
Since all the averages in both systems coincide, their free energies coincide, too. 
It is worth noting that the validity of the replacement of the operators ao and aj by c-numbers 
in the thermodynamic limit was confirmed in the numerous subsequent publications of various 
authors. Thus Bogoliubov's 1947 analysis of the many-body Hamiltonian by means of a c-number 
substitution for the most relevant operators in the problem, the zero-momentum mode operators, 
was justified rigorously. Since the Bogoliubov's 1947 analysis is one of the key developments in 
the theory of the Bose gas, especially the theory of the low density gases currently at the forefront 
of experiment H3HS1 this result is of importance for the legitimation of that theory. Additional 
arguments were given in study, where the Bose-Einstein condensation and spontaneous £7(1) sym- 
metry breaking were investigated on the basis of Bogoliubov's truncated Hamiltonian Hb for a 
weakly interacting Bose system, and adding a £7(1) symmetry breaking term y/V(\ao + A*Oq) 
to Hb, It was shown also, by using the coherent state theory and the mean- field approximation 
rather than the c-number approximations, that the Bose-Einstein condensation occurs if and only 
if the £7(1) symmetry of the system is spontaneously broken. The real ground state energy and 
the justification of the Bogoliubov c-number substitution were given by solving the Schroedinger 
eigenvalue equation and using the self-consistent condition. Thus the Bogoliubov c-number sub- 
stitutions were fully correct and the symmetry breaking causes the displacement of the condensate 
state. 

The concept of quasiaverages was introduced by Bogoliubov on the basis of an analysis of many- 
particle systems with a degenerate statistical equilibrium state. Such states are inherent to various 
physical many-particle systems Those are liquid helium in the superfluid phase, metals in the 
superconducting state, magnets in the ferromagnetically ordered state, liquid crystal states, the 
states of superfluid nuclear matter, etc. (for a review, see In case of superconductivity, 

the source v ^ fc v(k)(a\^a_ k < + a_^a^) was inserted in the BCS-Bogoliubov Hamiltonian, and 
the quasiaverages were defined by use of the Hamiltonian H v . In the general case, the sources are 
introduced to remove degeneracy. If infinitesimal sources give infinitely small contributions to the 
averages, then this means that the corresponding degeneracy is absent, and there is no reason to 
insert sources in the Hamiltonian. Otherwise, the degeneracy takes place, and it is removed by 
the sources. The ordinary averages can be obtained from quasiaverages by averaging with respect 
to the parameters that characterize the degeneracy. 
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N. N. Bogoliubov, Jr. 15 considered some features of quasiaverages for model systems with four- 
fermion interaction. He discussed the treatment of certain three-dimensional model systems which 
can be solved exactly. For this aim a new effective way of defining quasiaverages for the systems 
under consideration was proposed. 

Peletminskii and SokolovskiP^ have found general expressions for the operators of the flux den- 
sities of physical variables in terms of the density operators of these variables. The method of 
quasiaverages and the expressions found for the flux operators were used to obtain the averages 
of these operators in terms of the thermodynamic potential in a state of statistical equilibrium of 
a superfluid liquid. 

Vozyakov 4 ^ reformulated the theory of quantum crystals in terms of quasiaverages. He analyzed 
a Bose system with periodic distribution of particles which simulates an ensemble in which the 
particles cannot be regarded as vibrating independently about a position of equilibrium lattice 
sites. With allowance for macroscopic filling of the states corresponding to the distinguished sym- 
metry, a calculation was made of an excitation spectrum in which there exists a collective branch 
of gapless type. 

Peregoudov 4 ^ discussed the effective potential method, used in quantum field theory to study 
spontaneous symmetry breakdown, from the point of view of Bogoliubov's quasiaveraging proce- 
dure. It was shown that the effective potential method is a disguised type of this procedure. The 
catastrophe theory approach to the study of phase transitions was discussed and the existence 
of the potentials used in that approach was proved from the statistical point of view. It was 
shown that in the ease of broken symmetry, the nonconvex effective potential is not a Legendre 
transform of the generating functional for connected Green's functions. Instead, it is a part of 
the potential used in catastrophe theory. The relationship between the effective potential and 
the Legendre transform of the generating functional for connected Green's functions is given by 
Maxwell's rule. A rigorous rule for evaluating quasiaveraged quantities within the framework of 
the effective potential method was established. 

N. N. Bogoliubov, Jr. with M. Yu. Kovalevsky and co-authors^ developed a statistical approach 
for solving the problem of classification of equilibrium states in condensed media with sponta- 
neously broken symmetry based on the quasiaverage concept. Classification of equilibrium states 
of condensed media with spontaneously broken symmetry was carried out. The generators of 
residual and spatial symmetries were introduced and equations of classification for the order pa- 
rameter has been found. Conditions of residual symmetry and spatial symmetry were formulated. 
The connection between these symmetry conditions and equilibrium states of various media with 
tensor order parameter was found out. An analytical solution of the problem of classification of 
equilibrium states for superfluid media, liquid crystals and magnets with tensor order parameters 
was obtained. Superfluid 3 He, liquid crystals, quadrupolar magnetics were considered in detail. 
Possible homogeneous and heterogeneous states were found out. Discrete and continuous ther- 
modynamic parameters, which define an equilibrium state, allowable form of order parameter, 
residual symmetry, and spatial symmetry generators were established. This approach, which is 
alternative to the well-known Ginzburg-Landau method, does not contain any model assumptions 
concerning the form of the free energy as functional of the order parameter and does not employ 
the requirement of temperature closeness to the point of phase transition. For all investigated 
cases they found the structure of the order parameters and the explicit forms of generators of 
residual and spatial symmetries. Under the certain restrictions they established the form of the 
order parameters in case of spins 0, 1/2, 1 and proposed the physical interpretation of the studied 
degenerate states of condensed media. 
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3 Quantum Protectorate 



The "quantum protectorate" concept was formulated in paperP^l Its inventors, R. Laughlin and 
D. Pines, discussed the most fundamental principles of matter description in the widest sense of 
this word. They formulated their main thesis: emergent physical phenomena, which are regulated 
by higher physical principles, have a certain property, typical for these phenomena only. This 
property is their insensitivity to microscopic description. For instance, the crystalline state is 
the simplest known example of a quantum protectorate, a stable state of matter whose generic 
low-energy properties are determined by a higher organizing principle and nothing else. There are 
many other examples.^ These quantum protectorates, with their associated emergent behavior, 
provide us with explicit demonstrations that the underlying microscopic theory can easily have 
no measurable consequences whatsoever at low energies. The nature of the underlying the- 
ory is unknowable until one raises the energy scale sufficiently to escape protection. 
The existence of two scales, the low-energy and high-energy scales, relevant to the description 
of magnetic phenomena was stressed by the author of the present work in the papersj^E^l 
which were devoted to comparative analysis of localized and band models of quantum theory of 
magnetism. It was shown there, that the low-energy spectrum of magnetic excitations in the 
magnetically-ordered solid bodies corresponds to a hydrodynamic pole (k, cj — > 0) in the general- 
ized spin susceptibility x> which is present in the Heisenberg, Hubbard, and the combined s — d 
model. In the Stoner band model the hydrodynamic pole is absent, there are no spin waves there. 
At the same time, the Stoner single-particle's excitations are absent in the Heisenberg model's 
spectrum. The Hubbard model with narrow energy bands contains both types of excitations: the 
collective spin waves (the low-energy spectrum) and Stoner single-particle's excitations (the high- 
energy spectrum). This is a big advantage and flexibility of the Hubbard model in comparison 
to the Heisenberg model. The latter, nevertheless, is a very good approximation to the realistic 
behavior in the limit k,u —> 0, the domain where the hydrodynamic description is applicable, 
that is, for long wavelengths and low energies. The quantum protectorate concept was applied 
to the quantum theory of magnetism by the the present author in paperP where a criterion of 
applicability of models of the quantum theory of magnetism to description of concrete substances 
was formulated. The criterion is based on the analysis of the model's low-energy and high-energy 
spectra. There the detailed analysis was carried out of the idea of quantum protectorate^ in the 
context of quantum theory of magnetismPS It was suggested that the difficulties in the formu- 
lation of quantum theory of magnetism at the microscopic level, that are related to the choice 
of relevant models, can be understood better in the light of the quantum protectorate concept.^ 
We argued that the difficulties in the formulation of adequate microscopic models of electron and 
magnetic properties of materials are intimately related to dual, itinerant and localized behavior 
of electrons.^ We formulated a criterion of what basic picture describes best this dual behavior. 
The main suggestion was that quasiparticle excitation spectra might provide distinctive signatures 
and good criteria for the appropriate choice of the relevant model. A broad class of the problems 
of condensed matter physiciP^^ in the fields of the magnetism and superconductivity of complex 
materials were reconsidered in relation to these ideas. 



4 Irreducible Green Functions Method 

It was shown above that it becomes clear that only a thorough experimental and theoretical 
investigation of quasiparticle many-body dynamics of the many particle systems can provide the 
answer on the relevant microscopic picture. Many-particle systems where the interaction is strong 
have often complicated behavior, and require nonperturbative approaches to treat their properties. 
There are many different approaches to construction of generalized mean-field approximations; 
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however, all of them have a special-case character. The method of irreducible Green functions 
(IGF) allows one to tackle this problem in a more systematic fashion. In order to clarify this 
statement let us consider briefly the main ideas of the IGF approach that allows one to describe 
completely quasiparticle spectra with damping in a very natural way. When working with infinite 
hierarchies of equations for Green functions the main problem is finding the methods for their 
efficient decoupling, with the aim of obtaining a closed system of equations, which determine the 
Green functions. A decoupling approximation must be chosen individually for every particular 
problem, taking into account its character. This "individual approach" is the source of critique for 
being too not transparent, which sometimes appear in the papers using the causal Green functions 
and diagram technique. However, the ambiguities are also present in the diagram technique, when 
the choice of an appropriate approximation is made there. The decision, which diagrams one has 
to sum up, is obvious only for a narrow range of relatively simple problems. In the paper^ devoted 
to Bose-systems, and in the papers by the author of present work^EM^E] devoted to Fermi systems 
it was shown that for a wide range of problems in statistical mechanics and theory of condensed 
matter one can outline a fairly systematic recipe for constructing approximate solutions in the 
framework of irreducible Green's functions method. Within this approach one can look from a 
unified point of view at the main problems of fundamental characters arising in the method of 
two-time temperature Green functions. The method of irreducible Green functions is a useful 
reformulation of the ordinary Bogoliubov-Tyablikov method of equations of motionP^S 
We reformulated the two-time Green functions met ho j J^PoM UlMl ^ form which is especially 
adjusted to correlated fermion systems on a lattice and systems with complex spectra. A very 
important concept of the whole method is the generalized mean fields (GMFs), as it was formulated 
mJ22]|38] These GMFs have a complicated structure for a strongly correlated case and complex 
spectra, and are not reduced to the functional of mean densities of the electrons or spins when 
one calculates excitation spectra at finite temperatures. 

To clarify the foregoing, let us consider a retarded Green function of the formal 

G r = ((A(t),A\t'))) = -z8(t - t')([A(t)A^t%), V = ±1. (4.1) 

As an introduction to the concept of IGF, let us describe the main ideas of this approach in a 
symbolic and simplified form. To calculate the retarded Green function G(t — t'), let us write 
down the equation of motion for it 

uG(u) = ([A, A%) + (([A, | At)) w . (4.2) 

Here we use the notation ((A(t), A'(t'))) for the time-dependent Green function and ((A \ A^))^ 
for its Fourier transform.^ The notation [A, B]„ refers to commutation and anticommutation, 
depending on the value of r\ = ±. The essence of the method is as follows.^ It is based on the 
notion of the "IRRED U CIBLE" parts of Green functions (or the irreducible parts of the operators, 
A and A\ out of which the Green function is constructed) in terms of which it is possible, without 
recourse to a truncation of the hierarchy of equations for the Green functions, to write down the 
exact Dyson equation and to obtain an exact analytic representation for the self-energy operator. 
By definition, we introduce the irreducible part (ir) of the Green function 

W(([A,fl]_|At» = (([A,H]- - zA\A% (4.3) 

The unknown constant z is defined by the condition (or constraint) 

([[A,H]^\A%) = 0, (4.4) 

which is an analogue of the orthogonality condition in the Mori formalism EHEU us emphasize 
that due to the complete equivalence of the definition of the irreducible parts for the Green 
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functions { {ir) (([A, J5T]_|A+») and operators (^[A, = (L4, #]_) (ir) we will use both the 

notation freely ( is the same as <((A)( ir )|.B)) ). A choice one notation over another is 

determined by the brevity and clarity of notation only. From the condition (|4.4p one can find 

_ ([[A,H]^A%) _ My 

([A,A%) "Mo" (4 ' 5j 

Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, the 
irreducible Green functions are defined so that they cannot be reduced to the lower-order ones by 
any kind of decoupling. It is worth noting that the term "irreducible" in a group theory means a 
representation of a symmetry operation that cannot be expressed in terms of lower dimensional 
representations. Irreducible (or connected ) correlation functions are known in statistical mechan- 
ics. In the diagrammatic approach, the irreducible vertices are defined as graphs that do not 
contain inner parts connected by the G°-line. With the aid of the definition (14.3|) these concepts 
are expressed in terms of retarded and advanced Green functions. The procedure extracts all 
relevant (for the problem under consideration) mean-field contributions and puts them into the 
generalized mean-field Green function which is defined here as 

G0( W ) = y» At U ( 4. 6) 

To calculate the IGF (([A, B\-(t), A^t'))) in (|4.2p , we have to write the equation of motion for 
it after differentiation with respect to the second time variable t'. The condition of orthogonality 
(14. 4p removes the inhomogeneous term from this equation and is a very crucial point of the whole 
approach. If one introduces the irreducible part for the right-hand side operator, as discussed above 
for the "left" operator, the equation of motion (|4.2p can be exactly rewritten in the following form: 



G = G° + G°PG°. (4.7) 

The scattering operator P is given by 

P = (M )- 1 ( iir) (([A,HU[A\HU))^)(M )- 1 . (4.8) 

The structure of equation ( I4.8P enables us to determine the self-energy operator M by analogy 
with the diagram technique 

P = M + MG°P. (4.9) 

We used here the notation M for self-energy (mass operator in quantum field theory). From 
the definition (I4.9P it follows that the self-energy operator M is defined as a proper (in the 
diagrammatic language, "connected") part of the scattering operator M = (P) p - As a result, we 
obtain the exact Dyson equation for the thermodynamic double-time Green functions 

G = G° + G°MG. (4.10) 

The difference between P and M can be regarded as two different solutions of two integral equa- 
tions (|4,7p and (|4.10p . However, from the Dyson equation (|4.10p only the full GF is seen to be 
expressed formal solution of the form 

G= [(G )- 1 -M]-\ (4.11) 

Equation (|4.1ip can be regarded as an alternative form of the Dyson equation (|4.10p and the 
definition of M provides that the generalized mean-field GF G° is specified. On the contrary , for 
the scattering operator P, instead of the property + G M = 1, one has the property 

(G )- 1 -G- 1 = PG°G-\ 
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Thus, the very functional form of the formal solution (|4.1ip precisely determines the difference 
between P and M. 

Thus, by introducing irreducible parts of GF (or irreducible parts of the operators, out of which 
the GF is constructed) the equation of motion (|4.2p for the GF can exactly be (but using the 
orthogonality constraint (|4.4p ) transformed into the Dyson equation for the double-time thermal 
GF (I4.10p . This result is very remarkable because the traditional form of the GF method does 
not include this point. Notice that all quantities thus considered are exact. Approximations 
can be generated not by truncating the set of coupled equations of motions but by a specific 
approximation of the functional form of the mass operator M within a self-consistent scheme 
expressing M in terms of the initial GF 

M « F[G]. 

Different approximations are relevant to different physical situations. The projection operator 
technique has essentially the same philosophy. But with using the constraint ()4.4p in our ap- 
proach we emphasize the fundamental and central role of the Dyson equation for calculation of 
single-particle properties of many-body systems. The problem of reducing the whole hierarchy of 
equations involving higher-order GFs by a coupled nonlinear set of integro-differential equations 
connecting the single-particle GF to the self-energy operator is rather nontrivial. A characteristic 
feature of these equations is that besides the single-particle GF they involve also higher-order 
GF. The irreducible counterparts of the GFs, vertex functions, serve to identify correctly the 
self-energy as 

M = G 1 - G~ x . 

The integral form of the Dyson equation (|4.10p gives M the physical meaning of a nonlocal and 
energy-dependent effective single-particle potential. This meaning can be verified for the exact 
self-energy using the diagrammatic expansion for the causal GF. 

It is important to note that for the retarded and advanced GFs, the notion of the proper part 
M = (P) p is symbolic in natureP^EH j n a certain sense, it is possible to say that it is defined 
here by analogy with the irreducible many-particle T-matrix. Furthermore, by analogy with the 
diagrammatic technique, we can also introduce the proper part defined as a solution to the integral 
equation (|4.9p . These analogues allow us to better understand the formal structure of the Dyson 
equation for the double-time thermal GF, but only in a symbolic form . However, because of the 
identical form of the equations for GFs for all three types (advanced, retarded, and causal), we 
can convert our calculations to causal GF at each stage of calculations and, thereby, confirm the 
substantiated nature of definition (|4.9p . We therefore should speak of an analogy of the Dyson 
equation. Hereafter, we drop this stipulating, since it does not cause any misunderstanding. In a 
sense, the IGF method is a variant of the Gram-Schmidt orthogonalization procedure.^ 
It should be emphasized that the scheme presented above gives just a general idea of the IGF 
method. A more exact explanation why one should not introduce the approximation already in P, 
instead of having to work out M, is given below when working out the application of the method 
to specific problems. 

The general philosophy of the IGF method is in the separation and identification of elastic scat- 
tering effects and inelastic ones. This latter point is quite often underestimated, and both effects 
are mixed. However, as far as the right definition of quasiparticle damping is concerned, the 
separation of elastic and inelastic scattering processes is believed to be crucially important for 
many-body systems with complicated spectra and strong interaction. 

From a technical point of view, the elastic GMF renormalizations can exhibit quite a nontrivial 
structure. To obtain this structure correctly, one should construct the full GF from the complete 
algebra of relevant operators and develop a special projection procedure for higher-order GFs, in 
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accordance with a given algebra. Then a natural question arises how to select the relevant set of 
operators {Ai,^, -.-An} describing the "relevant degrees of freedom". The above consideration 
suggests an intuitive and heuristic way to the suitable procedure as arising from an infinite chain 
of equations of motion (14. 2p . Let us consider the column 

MA 

A 2 



\A n J 



where 



A 1 = A, A 2 = [A,H], A 3 = [[A,H],H},...A n = [[...[A,H}...H] 

n 

Then the most general possible Green function can be expressed as a matrix 



G = (( 



A 2 
\A n J 



a\ a\ 



At))) 



This generalized Green function describes the one-, two-, and n-particle dynamics. The equa- 
tion of motion for it includes, as a particular case, the Dyson equation for single-particle Green 
function, and the Bethe-Salpeter equation which is the equation of motion for the two-particle 
Green function and which is an analogue of the Dyson equation, etc . The corresponding reduced 
equations should be extracted from the equation of motion for the generalized Green function 
with the aid of special techniques such as the projection method and similar techniques. This 
must be a final goal towards a real understanding of the true many-body dynamics. At this 
point, it is worthwhile to underline that the above discussion is a heuristic scheme only, but not 
a straightforward recipe. The specific method of introducing the IGFs depends on the form of 
operators A n , the type of the Hamiltonian, and conditions of the problem. Here a sketchy form 
of the IGF method was presented. The aim was to introduce the general scheme and to lay the 
groundwork for generalizations. We demonstrated irp2ll2SI that the IGF method is a powerful tool 
for describing the quasiparticle excitation spectra, allowing a deeper understanding of elastic and 
inelastic quasiparticle scattering effects and the corresponding aspects of damping and finite life- 
times. In a certain sense, it provides a clear link between the equation-of-motion approach and the 
diagrammatic methods due to derivation of the Dyson equation. Moreover, due to the fact that it 
allows the approximate treatment of the self-energy effects on a final stage, it yields a systematic 
way of the construction of approximate solutions. 



5 Effective and Generalized Mean Fields 

The most common technique for studying the subject of interacting many-particle systems is 
to use the mean field theory. This approximation was especially popular in the theory of mag- 
netismP^^S To calculate the susceptibility and other characteristic functions of a system of 
localized magnetic moments, with a given interaction Hamiltonian, the approximation, termed 
the "molecular field approximation" was used widely. However, it is not an easy task to give the 
formal unified definition what the mean field is. In a sense, the mean field is the umbrella term 
for a variety of theoretical methods of reducing the many-particle problem to the single-particle 
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one. Mean field theory, that approximates the behavior of a system by ignoring the effect of fluc- 
tuations and those spin correlations which dominate the collective properties of the ferromagnet 
usually provides a starting and estimating point only, for studying phase transitions. The mean 
field theories miss important features of the dynamics of a system. The main intention of the 
mean field theories, starting from the works of J. D. van der Waals and P. Weiss, is to take into 
account the cooperative behavior of a large number of particles. It is well known that earlier 
theories of phase transitions based on the ideas of van der Waals and Weiss lead to predictions 
which are qualitatively at variance with results of measurements near the critical point. Other 
variants of simplified mean field theories such as the Hartree-Fock theory for electrons in atoms, 
lead to discrepancies of various kinds too. It is therefore natural to analyze the reasons for such 
drawbacks of earlier variants of the mean field theories. 

A number of effective field theories which are improved versions of the "molecular field approx- 
imation" were proposed. In our papers^K^^^ we stressed a specificity of strongly correlated 
many-particle systems on a lattice contrary to continuum (uniform) systems. The earlier concepts 
of molecular field were described in terms of a functional of mean magnetic moments (in magnetic 
terminology) or mean particle densities. The corresponding mean-field functional F[(n), {S z )] de- 
scribes the uniform mean field. Actually, the Weiss model was not based on discrete "spins" as 
is well known, but the uniformity of the mean internal field was the most essential feature of the 
model. In the modern language, one should assume that the interaction between atomic spins 
Si and its neighbors is equivalent to a mean (or molecular) field, M{ = Xo[h\ e ^ + hf" 1 ] and 
that the molecular field h\ m ^ is of the form h^ m ^ = J2i J{Rji)(Si) (above T c ). Here h ext is an 
applied conjugate field, xo is the response function, and J(Rji) is an interaction. In other words, 
the mean field approximation reduces the many-particle problem to a single-site problem in which 
a magnetic moment at any site can be either parallel or antiparallel to the total magnetic field 
composed of the applied field and the molecular field. The average interaction of % neighbors was 
taken into account only, and the fluctuations were neglected. One particular example, where the 
mean field theory works relatively well is the homogeneous structural phase transitions; in this case 
the fluctuations are confined in phase space. The next important step was made by L. Neel. He 
conjectured that the Weiss internal field might be either positive or negative in sign. In the latter 
case, he showed that below a critical temperature (Neel temperature) an ordered arrangement of 
equal numbers of oppositely directed atomic moments could be energetically favorable. This new 
magnetic structure was termed antiferromagnetism. It was conjectured that the two-sublattice 
Neel (classical) ground state is formed by local staggered internal mean fields. 
There is a number of the "correlated effective field" theories, that tend to repair the limitations 
of simplified mean field theories. The remarkable and ingenious one is the Onsager "reaction field 
approximation". He suggested that the part of the molecular field on a given dipole moment 
which comes from the reaction of neighboring molecules to the instantaneous orientation of the 
moment should not be included into the effective orienting field. This "reaction field" simply fol- 
lows the motion of the moment and thus does not favor one orientation over another (for details 
see RefePEHI). 

It is knowrJSSl that mean-field approximations, for example the molecular field approximation for 
a spin system, the Hartree-Fock approximation and the BCS-Bogoliubov approximation for an 
electron system are universally formulated by the Bogoliubov inequality: 

- p-HnfTret-W) < 

-p-MTveM* >) + Trei .; Hmf) l - (5-1) 

Here F is the free energy, and H m f is a "trial" or a "mean field" approximating Hamiltonian. 
This inequality gives the upper bound of the free energy of a many-body system. It is important 
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to emphasize that the BCS-Bogoliubov theory of superconductivity^"^^ was formulated on 
the basis of a trial Hamiltonian which consists of a quadratic form of creation and annihilation 
operators, including "anomalous" (off-diagonal) averages. The functional of the mean field (for 
the superconducting single-band Hubbard model) is of the following formP3 



1 M-A.) -<4^> 1 



The "anomalous" off-diagonal terms fix the relevant BCS-Bogoliubov vacuum and select the ap- 
propriate set of solutions. From the point of view of quantum many-body theory, the problem of 
adequate introduction of mean fields for system of many interacting particles can be most consis- 
tently investigated in the framework of the IGF method. A correct calculation of the quasiparticle 
spectra and their damping, particularly, for systems with a complicated spectrum and strong in- 
teraction^ reveals, that the generalized mean fields can have very complicated structure which 
cannot be described by a functional of the mean-particle density. 

To illustrate the actual distinction of description of the generalized mean field in the equation-of- 
motion method for the double-time Green functions, let us compare the two approaches, namely, 
that of Tyablikov-^ and of C alien P*^ We shall consider the Green function ((S + \S~)} for the 
isotropic Heisenberg model 

H = ~J2j(i-j)S i S j . (5.3) 

The equation of motion for the spin Green function is of the form 

u((S+\Sr)) u = (5.4) 
2{S z )5 l3 + Y J J(i-g){{StS z g - Spt\Sj)) w . 

9 

The Tyablikov decoupling expresses the second-order Green function in terms of the first (initial) 
Green function: 

((S+S z g \Sr)) = (S z )((S+\Sr)). (5.5) 
This approximation is an RPA-type; it does not lead to the damping of spin wave excitations 

E(q) = Y,J( i -9)(S z )ex P [i(R i -R g )q\=2(S z )(Jo-J q ). (5.6) 

9 

The reason for this is rather transparent. This decoupling does not take into account the inelastic 
magnon-magnon scattering processes. In a sense, the Tyablikov approximation consists of approx- 
imating the commutation relations of spin operators to the extent of replacing the commutation 
relation = 2S?5 tJ by [S+,Sr}_ = 2(S z )5 ij . 

CallerPSl has proposed an improved decoupling approximation in the method of Tyablikov in the 
following form: 

(S z )((SpB)) - a(S g Sp((S+\B)). 

Here < a < 1. To clarify this point, it should be reminded that for spin 1/2 ( the procedure was 
generalized by Callen to an arbitrary spin), the spin operator S z can be written as Sg = S — S~S+ 
or Sg = ^(S^S~ — 5 ~ 5+ ) . It is easy to show that 

Sg =aS H — S g S g — Sg ^ . 
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The operator S g S£ represents the deviation of (S z ) from S. In the low-temperature region, this 
deviation is small, and a ~ 1. Similarly, the operator ^(S^S~ — S~S£) represents the deviation 
of (S z ) from 0. Thus, when (S z ) approaches to zero, one can expect that a ~ 0. Thus, in this 
way, it is possible to obtain a correction to the Tyablikov decoupling with either a positive or 
negative sign, or no correction at all, or any intermediate value, depending on the choice of a. 
The above Callen arguments are not rigorous , for, although the difference in the operators S + S~ 
and S~S + is small if (S z ) ~ 0, each operator makes a contribution of the order of S, and it is 
each operator which is treated approximately, not the difference. There are some other drawbacks 
of the Callen decoupling scheme. Nevertheless, the Callen decoupling was the first conceptual 
attempt to introduce the interpolation decoupling procedure. Let us note that the choice of a = 
over the entire temperature range is just the Tyablikov decoupling (|5.5p . 
The energy spectrum for the Callen decoupling is given by 

E(q) = 2(S Z )((J - J q ) + ^£[J(fc) - J{k - q)]N{E(k))). (5.7) 

k 

Here N{E{k)) is the Bose distribution function N(E(k)) = [exp(E(k)f3) -l] -1 . This is an implicit 
equation for N(E(k)), involving the unknown quantity {S z ) . For the latter an additional equation 
is givenPSl Thus, both these equations constitute a set of coupled equations which must be solved 
self-consistently for (S z ). This formulation of the Callen decoupling scheme displays explicitly 
the tendency of the improved description of the mean field. In a sense, it is possible to say 
that the Callen work dates really the idea of the generalized mean field within the equation-of- 
motion method for double-time GFs, however, in a semi-intuitive form. The next essential steps 
were made by PlakidcP^ for the Heisenberg ferromagnet and by Kuzemsky^ for the Hubbard 
model. Later many approximate schemes for decoupling the hierarchy of equations for GF were 
proposed, improving the Tyablikov's and Callen's decouplings. Various approaches generalizing 
the random phase's approximation in the ferromagnetism theory for wide ranges of temperature 
were considered in paper by Czachor and HolasP^ As was mentioned above, the correct definition 
of generalized mean fields depends on the condition of the problem, the strength of interaction, the 
choice of relevant operators, and on the symmetry requirements. The most important conclusion 
to be drawn from the present consideration is that the GMF, in principle, can have quite a 
nontrivial structure and cannot be reduced to the mean-density functional only. 

6 Quasiaverages and Irreducible Green Functions Method 

In condensed matter physics, the symmetry is important in classifying different phases and under- 
standing the phase transitions between them. There is an important distinction between the case 
where the broken symmetry is continuous (e.g. translation, rotation, gauge invariance) or discrete 
(e.g. inversion, time reversal symmetry) P2 The Goldstone theorem states that when a continuous 
symmetry is spontaneously broken and the interactions are short ranged a collective mode (excita- 
tion) exists with a gapless energy spectrum (i.e. the energy dispersion curve starts at zero energy 
and is continuous). Acoustical phonons in a crystal are prime examples of such so-called gapless 
Goldstone modes. Other examples are the Bogoliubov sound modes in (charge neutral) Bose 
condensates32HSl a nd spin waves (magnons) in ferro- and antiferromagnets. N. N. Bogoliubov and 
then Y. Nambu in in their works shown that the general features of superconductivity are in fact 
model independent consequences of the spontaneous breakdown of electromagnetic gauge invari- 
ance. It is important to emphasize that the BCS-Bogoliubov theory of superconductivitj^^IES] 
was formulated on the basis of a trial Hamiltonian which consists of a quadratic form of creation 
and annihilation operators, including "anomalous" (off-diagonal) averages.^ 3 The strong-coupling 
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BCS-Bogoliubov theory of superconductivity was formulated for the Hubbard model in the local- 
ized Wannier representation in RefsF^®"^ Therefore, instead of the algebra of the normal state's 
operator a^, a\ a and rii a , for description of superconducting states, one has to use a more general 
algebra, which includes the operators di a , a\ a , rii a and ai a ai- a , a\ a a\_ a . The relevant generalized 
one-electron Green function will have the following formPSl^^Ql 

(Gu G 12 \ ( {{aia\a) a )) {{a ia \aj- a ))\ 

As it was discussed in Refs.pM^the off-diagonal (anomalous) entries of the above matrix select 
the vacuum state of the system in the BCS-Bogoliubov form, and they are responsible for the 
presence of anomalous averages. For treating the problem we follow the general scheme of the 
irreducible Green functions methodJ22E^ In this approach we start from the equation of motion 
for the Green function Gij(ui) (normal and anomalous components) 

y^juSjj - tij)((a ja \al a )) = 5u, + (6.2) 
j 

U {{a i(T ni- a \al /a . crUn\ala))i 

nj 

X>^ + %)((a]->L>>= ( 6 - 3 ) 
i 

- U {{ a \-a n i°\ a \>o)) + ^2 V M(a]_ a u n \a], a }}. 

nj 

The irreducible Green functions are introduced by definition 

(^((oiaot^oi-alot^)),,,) = {{a ia a\_ a ai^ a \a\, a )) ^ - (6.4) 
-(ni- a )Gn + (aj (T a i _ (7 )((aJ_ (T |a t T , (T )) w , 

( M ((4* a i«4~*\4*))u) = ((4x a ^4->L>>" - 
-{n ia )G 2 i + (ata|_ (J )((a i(T |a l t , (T )) tJ . 

The self-consistent system of superconductivity equations follows from the Dyson equatiorP^^EB 

G tl ,(u) = G% (to) +Y J G Q ij {u)M jf {u)G j , i ,{u>). (6.5) 
jj' 

The mass operator Mjj>(uj) describes the processes of inelastic electron scattering on lattice vibra- 
tions. The elastic processes are described by the quantity (|5.2|) . Thus the "anomalous" off-diagonal 
terms fix the relevant BCS-Bogoliubov vacuum and select the appropriate set of solutions. The 
functional of the generalized mean field for the superconducting single-band Hubbard model is of 
the form Yf a . A remark about the BCS-Bogoliubov mean-field approach is instructive. Speaking 
in physical terms, this theory involves a condensation correctly, in spite that such a condensation 
cannot be obtained by an expansion in the effective interaction between electrons. Other mean 
field theories, e.g. the Weiss molecular field theory and the van der Waals theory of the liquid-gas 
transition are much less reliable. The reason why a mean-field theory of the superconductivity 
in the BCS-Bogoliubov form is successful would appear to be that the main correlations in metal 
are governed by the extreme degeneracy of the electron gas. The correlations due to the pair 
condensation, although they have dramatic effects, are weak (at least in the ordinary supercon- 
ductors) in comparison with the typical electron energies, and may be treated in an average way 
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with a reasonable accuracy. It should be emphasized that the high-temperature superconductors, 
discovered two decades ago, motivated an intensification of research in superconductivity, not only 
because applications are promising, but because they also represent a new state of matter that 
breaks certain fundamental symmetries. These are the broken symmetries of gauge (superconduc- 
tivity), reflection (d-wave superconducting order parameter), and time-reversal (ferromagnetism) . 
Superconductivity and antiferromagnetism are both the spontaneously broken symmetries. The 
question of symmetry breaking within the localized and band models of antiferromagnets was 
studied by the author of this work in Refs! 22 ^ 32 ^ 3 It has been found there that the concept of 
spontaneous symmetry breaking in the band model of magnetism is much more complicated than 
in the localized model. In the framework of the band model of magnetism one has to additionally 
consider the so called anomalous propagators of the form 

FM:G /m ~((a to ;4j), 
AFM : G afm ~ ((a k+Q(T ;al +Q , a ,)). 

In the case of the band antiferromagnet the ground state of the system corresponds to a spin- 
density wave (SDW), where a particle scattered on the internal inhomogeneous periodic field gains 
the momentum Q — Q' and changes its spin: a — >• & ' . The long-range order parameters are defined 
as follows: 

FM : m = l/N ^ a ka a h-a), (6.6) 

ka 

AFM : M Q = Y,(4* a k+Q-*)- (6.7) 

ka 

It is important to stress, that the long-range order parameters here are functionals of the internal 
field, which in turn is a function of the order parameter. Thus, in the cases of rotation and trans- 
lation invariant Hamiltonians of band ferro- and antiferromagnetics one has to add the following 
infinitesimal sources removing the degeneracy: 

FM : ufi B H x a L a k-a, (6.8) 
ka 

AFM : ufi B H 4a a k+Q-a- (6.9) 
kQ 

Here, v — )■ after the usual in statistical mechanics infinite-volume limit V — > oo. The ground 
state in the form of a spin-density wave was obtained for the first time by Overhauser. There, 
the vector Q is a measure of inhomogeneity or translation symmetry breaking in the system. The 
analysis performed by various authors showed that the antiferromagnetic and more complicated 
states (for instance, ferrimagnetic) can be described in the framework of a generalized mean-field 
approximation.^ In doing that we have to take into account both the normal averages {a\ a ai a ) 
and the anomalous averages (a\ a ai- a ). It is clear that the anomalous terms break the original 
rotational symmetry of the Hubbard Hamiltonian. Thus, the generalized mean-field's approxi- 
mation for the antiferromagnet has the following formp3 m- a ai a ~ (nj_ (7 )aj - — (a|_ cr aj -)aj_ cr . A 
self-consistent theory of band antiferromagnetisrrPH was developed by the author of this work us- 
ing the method of the irreducible Green functions ! 22 H 38 I The following definition of the irreducible 
Green functions was used: 

{{ a k+pa a \>J r q-a a q-a\ a \ a ))u] = {{ a k+pa a \,+q-a a q-a\ a \ u )) 10 ~ 

5 P fi(n q ~a)G kr7 - {a k+pa a\ >+q _ a ){{a q -- a \a\ (T )) w . (6.10) 
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The algebra of relevant operators must be chosen as follows ((cw, aj CT , n^, a\ a ai- a ). The corre- 
sponding initial GF will have the following matrix structure 



Qafm 



\a) a )) ((a la \a]_ a )) 



J) {{«-*\a}_J) 



The off-diagonal terms select the vacuum state of the band's antiferromagnet in the form of a spin- 
density wave. With this definition, one introduces the so-called anomalous (off-diagonal) Green 
functions which fix the relevant vacuum and select the proper symmetry broken solutions. The 
theory of the itinerant antiferromagnetisirpH was formulated by using sophisticated arguments of 
the irreducible Green functions method in complete analogy with our description of the Heisenberg 
antiferromagnet at finite temperatures.^ For the two-sublattice antiferromagnet we used the 
matrix Green function of the form 

G(k ^-{((s+\s Zka )) ((s+\s Zkb )))- t 6 - 11 ) 

Here, the Green functions on the main diagonal are the usual or normal Green functions, while the 
off-diagonal Green functions describe contributions from the so-called anomalous terms, analogous 
to the anomalous terms in the BCS-Bogoliubov superconductivity theory. The anomalous (or off- 
diagonal) average values in this case select the vacuum state of the system precisely in the form 
of the two-sublattice Neel stated 



7 Conclusions 

In the present work we shown that the development and improvement of the methods of quantum 
statistical mechanics still remains quite an important direction of research E3231 i n particular, the 
Bogoliubov's method of quasiaverages gives the deep foundation and clarification of the concept of 
broken symmetry.^ It makes the emphasis on the notion of a degeneracy and plays an important 
role in equilibrium statistical mechanics of many-particle systems. According to that concept, 
infinitely small perturbations can trigger macroscopic responses in the system if they break some 
symmetry and remove the related degeneracy (or quasidegeneracy) of the equilibrium state. As 
a result, they can produce macroscopic effects even when the perturbation magnitude is tend to 
zero, provided that happens after passing to the thermodynamic limit. 

We have discussed the theory of the correlation effects for many-particle interacting systems 
using the ideas of quasiaverages for interacting electron and spin systems on a lattice. The 
workable and self-consistent irreducible Green functions approach to the decoupling problem for 
the equation-of-motion method for double-time temperature Green functions has been presented. 
The main advantage of the formalism consists in the clear separation of the elastic scattering 
corrections (generalized mean fields) and inelastic scattering effects (damping and finite lifetimes). 
These effects could be self-consistently incorporated in a general and compact manner. Using 
the IGF method, it is possible to obtain a closed self-consistent set of equations determining the 
relevant Green functions and self-energy. These equations give a general microscopic description of 
correlation effects. Moreover, this approach gives the workable scheme for the definition of relevant 
generalized mean fields written in terms of appropriate correlators. This picture of interacting 
many-particle systems on a lattice is far richer and gives more possibilities for the analysis of 
phenomena which can actually take place. In this sense the approach we described produces more 
advanced physical picture of the quasiparticle many-body dynamics. 
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